Preprint typeset in JHEP style - HYPER VERSION 



Preprint 



Bulk Stabilization, the Extra-Dimensional 

Higgs Portal and Missing Energy in Higgs Events 



Ross Diener and CP. Burgess 

Department of Physics & Astronomy, McMaster University 
1280 Main Street West, Hamilton ON, Canada. 

Perimeter Institute for Theoretical Physics 

31 Caroline Street North, Waterloo ON, Canada. 



Abstract: To solve the hierarchy problem, extra-dimensional models must explain why the 
new dimensions stabilize to the right size, and the known mechanisms for doing so require bulk 
scalars that couple to the branes. Because of these couplings the energetics of dimensional 
stabilization competes with the energetics of the Higgs vacuum, with potentially observable 
effects. These effects are particularly strong for one or two extra dimensions because the bulk- 
Higgs couplings can then be super-renormalizable or dimensionless. Experimental reach for 
such extra-dimensional Higgs 'portals' are stronger than for gravitational couplings because 
they are less suppressed at low-energies. We compute how Higgs-bulk coupling through such 
a portal with two extra dimensions back-reacts onto properties of the Higgs boson. When the 
KK mass is smaller than the Higgs mass, mixing with KK modes results in an invisible Higgs 
decay width, missing-energy signals at high-energy colliders, and new mechanisms of energy 
loss in stars and supernovae. Astrophysical bounds turn out to be complementary to collider 
measurements, with observable LHC signals allowed by existing constraints. We comment 
on the changes to the Higgs mass-coupling relationship caused by Higgs-bulk mixing, and 
how the resulting modifications to the running of Higgs couplings alter vacuum-stability and 
triviality bounds. 
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1. Introduction 



In particle physics it is the best of times, and it is the worst of times. On one hand the recent 
discovery Q of a new particle at the LHC moves us into the long-awaited study of the new 
particle's properties, after several decades spent exploring the physics of constraints. If the 
new particle's interpretation as a Higgs — or the Higgs, if the Standard Model description 
continues to work — survives, then we can anticipate an unprecedented new era probing 
vacuum physics. 

On the other hand, the LHC has yet to produce compelling evidence for the kinds of 
physics widely expected to lie beyond the Standard Model. The hierarchy problem lies at 
the heart of these expectations, leading broadly to three main options^ for LHC-observable 
new electroweak physics over the years: compositeness models ||3|, supersymmetry ^ 
(linearly realized^); and extra-dimensional scenarios (both warped |^ and unwarped |llO|| ). 
Absent compelling evidence for any of these three categories, it is crucial for theorists to seek 
new ways to distinguish the mechanisms underlying each. 

The purpose of this paper is to identify new ways to use the properties of the Higgs 



to explore extra-dimensional models. Building on earlier work — in 5D |11] and higher 
dimensional scenarios |l3|, U, |l5| — we track how the vacuum energetics of the Higgs 
potential interacts with the physics that stabilizes the extra dimensions, and show how this 
can open a new observable portal onto extra-dimensional dynamics. 

At present most bounds on extra dimensions come from the kinematics of mixing and 
energy loss with the bulk gravitational degrees of freedom [^]. Yet a central part of solving 
the hierarchy problem using extra dimensions is understanding the vacuum physics that 
stabilizes their size at the required value, both for RS models (where the hierarchy between 
the electroweak and Planck scales comes from a size-dependent warp factor) and for ADD- 
type models (where it is the large extra-dimensional volume itself that provides the hierarchy) . 
All of the known mechanisms for this stabilization involve introducing new bulk degrees of 
freedom (typically scalar fields), whose couplings to ordinary matter are only slightly less 
robust than those of the metric. It is these couplings to the Higgs that we aim to constrain. 
We identify two kinds of observable consequences for these couplings. 

• Modified Higgs mass- coupling relations: due to the dependence of the Higgs potential 
on the new bulk fields. The interplay between these two fields changes the relationship 
between the Higgs mass and its couplings relative to Standard Model expectations; 

• Contributions to the Higgs invisible 'width': due to mixing between the Higgs and bulk 
states. In particular, we find that the expected LHC bounds on this width are compet- 
itive with bounds from lower-energy observables, such as energy loss from astrophysical 
systems, anomalous magnetic moments and the like. 



^These need not be mutually exclusive, with some composite models potentially being equivalent to some 
extra-dimensional models |^]. 

^See, however, for how supersymmetry could be present (but nonlinearly realized at electroweak 
energies and below, without requiring the existence of the superpartners that remain missing from experiments. 
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Two things are crucial about both of these effects. First, because the bulk fields involved 
are not the graviton, their couplings need not be precisely gravitational in strength. In 
particular (a point made earlier for ADD models in |jl^) depending on the number of extra 
dimensions present, they can involve dimensionless couplings, and so be less suppressed at 
low energies than are graviton interactions. (Dimensionless couplings can also arise for Higgs- 
curvature interactions in 6D, but unlike the Higgs-bulk portal they remain suppressed at low 
energies because of the derivative nature of the curvature couplings |l2|.) 

Second, the interplay between the (brane-localized) Higgs and extra-dimensional (bulk) 
stabilization mechanisms depends crucially on understanding how branes back-react on the 
bulk. Although this is understood relatively well for branes with one transverse dimension 



(such as arise in RS models) in terms of Israel junction conditions |17|, it has only recently 
been systematically developed for branes with two or more transverse dimensions, such 
as appear in the ADD picture. The understanding of codimension-2 back-reaction came 
comparatively late because of technical complications associated with the divergence of bulk 
fields near brane positions (which happens only with two or more transverse dimensions), and 
the need to absorb these into renormalizations of the brane couplings . 

1.1 Higher-dimensional stabilization 

Until recently a big competitive advantage of RS models over ADD models was the existence 
of a simple and robust way to stabilize the extra dimensions: the Goldberger-Wise mechanism 
[^ ]. In this mechanism a bulk scalar field is introduced that couples to the branes situated 
at both ends of the RS scenario's one extra dimension, with couplings chosen to frustrate the 
scalar's ability to reach a constant vacuum configuration. (This can be achieved by having 
branes disagree with one another about the field value that minimizes the scalar potential.) 
Because branes are located at specific places in the extra dimension, the resulting frustration 
sets up gradients in the bulk scalar that make the minimum energy depend on the distance 
between the branes (and so also on the extra-dimensional size). An attractive feature of 
the RS model is that the warp factor then naturally exponentiates a modestly large extra- 
dimensional size into an enormous electroweak hierarchy. (Similar frustration can also be 
arranged with bulk scalars in more than one extra dimension, with sometimes intriguing 



implications for the Higgs vacuum [ 13 , |T 



A similarly robust mechanism for stabilizing large dimensions has been missing for stan- 



dard ADD models, but an analogue was recently found |22] for their supersymmetric general- 



izations ||2^, 24, |25|] by applying to them a 6D cousin [^6| of the Goldberger-Wise mechanism. 
In such theories the extra dimensions are stabilized classically through flux-stabilization, as 
is often possible for supersymmetric systems (and for which 6D systems provided the first ex- 
amples [^4|). In this mechanism the flux of a bulk magnetic field (which is typically required 
by anomaly cancellation to exist among the field content of the 6D supergravity |^, p5[ ) 
threads the two extra dimensions, that have the topology of a sphere. Dirac quantization of 
this flux makes it energetically costly to shrink the dimensions, providing a counterbalance 
against its gravitational collapse. 
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However complete stabilization purely within the bulk is never quite possible because 
of a classical scale invariance of the 6D supergravity action, which leaves a flat direction 
parameterized by a bulk scalar field, x (the 'dilaton', which sits within the 'extended' metric 
supermultiplet). Flux stabilization relates the extra-dimensional radius to this flat direction 
through the expression 

where £ is a length of order (but, in controlled calculations, parametrically moderately larger 
than) the 6D Planck scale, set by the flux stabilization. 

Fixing r completely requires breaking the classical scale invariance, and lifting the clas- 
sical bulk flat direction. As shown in ref. |22|, this can be achieved classically through its 
couplings to branes, whose interactions need not share the scale invariance of the bulk. In 
particular it is not difficult to arrange for moderately large negative values. Once this is done 



flux stabilization — via eq. (1.1) — ensures the resulting radius is exponentially large in Xi 



naturally ensuring an exponentially large hierarchy in these models as well. In the special 



case of super symmetric ADD models p5| , 27| £ ~ (lOTeV)"^, and so micron-sized dimensions 
can be achieved with x ~ —70. But the stabilization mechanism itself doesn't rely on using 
an ADD framework, and could equally well apply if it were the Kaluza-Klein (KK) scale that 
were of electroweak size. 



Of course quantum effects can modify eq. (1.1) because they break the classical bulk scale 



invariance. But since each loop breaks scale invariance by a specific amount, these turn out 
to generate corrections as a series in e^-^ |^8|, and so do not ruin the exponentially large size 
of r. 

1.2 Relevance to the Higgs 

From the point of view of the Higgs, what is important about the above mechanisms (in both 
5 and 6 dimensions) is that they require the presence of a coupling between a bulk scalar 
field and the brane on which the Higgs sits. For instance, in the 6D case the most general 
renormalizable interactions between a brane-localized Standard Model and a (canonically 
normalized) electroweak singlet bulk scalar, <I>, have the form 



'int 



d^x U [h^H, , (1.2) 



with 



where 



= T + fil<^h + Y'^l - it'H - H^H + x(^H^Hy , (1.3) 

r:=ro + ^ + ^, fil:=X2V^ and fil := Xv\ (1.4) 
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and = ^{x, y = yb), denotes the evaluation of the bulk scalar at the position of the brane. 
It is the dimensionless coupling g that represents the unique Standard Model portal into extra 
dimensions within this six-dimensional context. 

This means that the vacuum energetics of the Higgs field interacts with the physics 
that stabilizes the extra dimensions, and both Higgs and bulk fields must be varied to find 
the proper vacuum configuration. In particular, the bulk scalar couplings can act to help 
or hinder the propensity for electroweak symmetry breaking. For instance, to the extent 
that large volume requires <I>fc < we see that this acts to increase^ the effective value 



f^'neS ~ A ~ 9^b-, and so assists the formation of a nonzero v.e.v. for H. In what follows § 2.2 
fleshes this out more explicitly, with care being taken to handle properly the renormalizations 
required because actually diverges at the brane position. 
Similarly, using the replacement 



H = ^\ J, (1-5) 




in the term gH^H^b contributes to Higgs-bulk mixing, and so to invisible channels where 
energy leaks into the extra dimensions during Higgs-production processes. As we show be- 
low, such leakage looks like a Higgs invisible width, and so is subject to similar constraints. 
Furthermore, since the coupling g is dimensionless, this loss rate is less suppressed at lower 
energies than would have been true for gravitational energy loss, and so allows better bounds 
and opportunities for detection §2.2 computes this more carefully, extending the re- 



sults of |15] by taking full account of the Higgs-KK mixing brought about by brane-bulk 
back-reaction. 

The calculation in 6D in many ways resembles earlier work which considered Higgs- 



curvature mixing |12|, of the form H^H R, but with three differences. First, because the cur- 
vature couplings involve more derivatives than do the Higgs-scalar couplings, the curvature 
mixing remains suppressed at low energies (like other gravitational interactions). Secondly, 
unlike these earlier calculations, we are able to compute both the real and imaginary parts of 
the Higgs production amplitude and so can compute the full line-shape rather than just its 
effective width. We can do so because our treatment of back-reaction allows us to renormalize 
the divergences that complicate obtaining the real part, associated with the near-brane diver- 
gences of the bulk fields. This technology allows us to extend the study of mixing to invisible 
final states in astrophysics, and at colliders. Finally, we include all possible renormalizable 
interactions, including in particular the quadratic self-coupling, A2, for the bulk field on the 
brane. This inclusion has important consequences, since bounds on g weaken with increasing 
A2, ultimately allowing a detectable invisible width at the LHC be consistent with strong 
constraints from low-energy astrophysics (see Fig. |l|). 

Although our results apply both to the cases of large dimensions {nih ^ m^K) and 
small ones (m/i <^ rnj^)^ when discussing the phenomenology we focus on the case when 

^For g > 0. 
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Figure 1: A plot summarizing the various constraints and discovery potential of Higgs-bulk mixing 
in the G — gj^fa vs A2 — X^jd plane (renormalized at f = l/rrih), in the large-volume limit {rrih ^ 
TOka')- The quantity a '-^ 1 is a measure of the defect angle near the brane, as defined in detail in 



2.1. The dark (blue) shaded region is the region disfavoured by LHC global fits. The medium (blue) 



shaded region is the conservative bound from nucleon-bulk bremsstrahlung in SN1987a, assuming 
Tgjv = 20 MeV. The lightest (gray) shade denotes regions excluded by demanding no Landau poles 
below /i* = 1 TeV, with the vertical dotted lines denoting how this bound changes with the choice of 
ultraviolet scale . Also plotted are lines of constant invisible branching ratio B that will be probed 
with additional data at the LHC or future experiments, all of which constrain this quantity. 



the dimensions are large. We find, as did earlier authors [12, |T^, that a Higgs undergoing 
Higgs-bulk mixing in many ways resembles a Higgs that can decay into invisible channels. 
Indeed once both real and imaginary parts of the amplitude are computed, we find that 
the resemblance becomes perfect for processes with the Higgs resonantly produced in the 
narrow- width limit. 

However, because the resonant, narrow-width limit is not always sufficient, there are also 
important differences between a bulk-mixed Higgs and one with access to invisible decays. 
Most important among these is the existence of strong bounds from astrophysical processes 
like SN1987a. These are not normally relevant for a Higgs with invisible decay channels (or 
for Higgs-bulk mixing through the H R term), because the rate for producing the Higgs 
is too small at low energies to give an appreciable energy-loss channel. The same is not 
true for Higgs-bulk mixing in the scalar potential, however, since this is not suppressed at 
low energies, and is not dominated by resonant Higgs production. It is instead enhanced 
by the kinematic availability of a large number of very light states for which the couplings 
cannot be neglected. The resulting constraint is shown in Fig. 0, together with the constraint 
coming from the successful Standard Model description of the observed Higgs, and contours 
indicating the size of the effective invisible Higgs width. Although astrophysics furnishes a 
very strong constraint, it does not exclude the range of interest to future LHC measurements. 
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It does not do so because it is not a resonant process, and so involves a different combination 
of parameters than are measured at the LHC in the g — X2 plane. 

The next sections present the details of this analysis as follows. First, §|2| derives the 
main expressions for the propagation eigenstates when the Higgs mixes with bulk. Armed 
with these calculations we discuss some of the resulting phenomenology in §|3|, including the 
lineshape for Higgs production at the LHC, and various bounds from lower energy phenomena. 
Our conclusions are summarized briefly in §13. 



2. Higgs-Bulk Dynamics 

In this section we compute in detail the implications of a Higgs-bulk interaction of the form 
given in eq. For concreteness we restrict from here on to the 6D case, for which back- 

reaction issues are much less well-explored. 



We start, in §2.1 where the bulk and brane actions and field equations, including the 



conditions for back-reaction, are described. The next subsection, §2.2 then calculates how 
eq. ( |1 . 1| ) changes the energy minimization for the Higgs and bulk scalar fields and so alters the 
expression for the Higgs v.e.v. in terms of the parameters in its potential. This is followed 
in §2.3 by a calculation of the spectrum of fluctuations, including a treatment of how the 
on-brane Higgs mixes with the bulk KK states. We consider two limits of interest in this 
mixing, depending on whether the KK mass is small or of the same order as the on-brane 
Higgs mass. The former is of most interest for ADD and super symmetric large-dimension 
(SLED) scenarios, while the latter would be of interest for dimensions whose KK scale is of 
order the electroweak scale. We specialize to the case of large dimensions when examining 
phenomenology more explicitly. 



2.1 Field equations and back- reaction 

We start by describing the 6D bulk and 4D brane systems of interest. For simplicity we 
focus purely on a single bulk scalar field, coupled to a Standard Model Higgs doublet on a 
space-filling codimension-2 brane situated at a specific spot in the two extra dimensions. 



The Action 

Consider a massless, free, 6D bulk scalar, with action 

Sb = -J d'^xV^ Q G^'^'d,,^ 9^$) . (2.1) 

We do not include a scalar potential in the bulk, and for ADD-type models this could be natu- 
rally enforced through a shift symmetry. The presence of scalars in the gravity supermultiplet 
and in the massless hypermultiplet representations of 6D supersymmetry also make it natural 
include light bulk scalars when the extra dimensions are supersymmetric. In the simplest case 
for bulk stabilization the supergravity of interest is gauged, chiral supergravity, and <I> = V'^x 
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represents the canonically normalized dilaton that transforms in the (extended) gravity mul- 
tiplet In this case there is a bulk scalar potential, Ub{x) oc e^, which considerably 

complicates the treatment of fluctuations once the metric is included. However because the 
gravitational couplings are RG-irrelevant we omit them for simplicity of presentation, and 
expect our considerations explored here to apply at sufficiently low energies. 

Next, consider a space- filling 4D brane that is located at a particular point, y = yb, within 
the extra dimensions. With eq. ( |1.3| ) in mind we take the brane action to be 

Sb = I d^xV^ (^Cs,, - To - ^ (cDb + V^f - g H^H , (2.2) 

where $b := $(x,y = yh) and 7^1, = GMN{x,y = yb)d^z'^di,z'^ is the induced metric on the 
brane, whose world-sheet is denoted = {x'^,?/'" = y^{x)}. Csm denotes the Standard 
Model action, but for the present purposes we need only work with its Higgs part: 

(2 \ 2 
H^H-i^] . (2.3) 

Thus, the complete, on-brane scalar potential reads 

Ub = T- fxlH^H + X{H^H)^ + fil<^b + Y'^l+9 H <^b , (2.4) 

as anticipated in eq. (^]^). This contains all possible terms involving only H and the Standard 
Model that are local and involve only relevant or marginal couplings. 

Background Geometry 

For the purposes of discussing Higgs energetics, consider the following unwarped, axisymmet- 
ric background geometry, 

ds^ = g,,„ dx*Mx^ = r/^^ dx^'dx" + f{r) dO^ + dr^ , (2.5) 

where r denotes proper distance away from the brane on which the Higgs resides. We allow 
for the possibility of a conical singularity at this brane by allowing a defect angle: < 9 < 
27ra, with < q < 1. Control of approximations usually requires a small defect angle, so 
|q — 1| ^ 1. In real examples of interest the radial coordinate runs through a finite range, 
< r < ttR, with r = ttR associated with another 4D brane at the opposite end of the extra 
dimensions. 

For the present purposes we ask for simplicity that the singular behaviour of the extra- 
dimensional geometry be no worse than a conical singularity at the brane position, and so 
require /(r) ~ r for r <^ R. This is not the most general case but is broad enough to include 
a variety of back-reacted examples, such as locally flat extra dimensions — corresponding to 
/(r) = r — and spherical (or rugby-ball, for nonzero deficit angle) extra dimensions — for 
which /(r) = Rsm{r /R) — as well as other potentially more exotic geometries. 
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We do not specify /(r) explicitly other than this near-brane limit. This generality is 
possible because for collider applications to ADD- type models not much depends on f{r). 
Physically, this is because it is only the enormous phase space associated with the large 
number of very high energy modes that allows observably large contributions to collider 
physics at all. But these modes have such short wavelengths that they are insensitive to the 



large-scale shape of the extra dimensions (see, for example, |29] for explicit calculations that 
illustrate this point). 

2.2 Vacuum configurations 

We now seek vacuum solutions to the coupled brane-bulk field equations, subject to the 
assumptions of 4D Lorentz invariance and axisymmetry in the extra dimensions. 

Bulk field equations and vacuum solutions 

Using $ = $(r) in the bulk field scalar equation, = 0, then gives 

dr{fdr<^) = 0, (2.6) 

which integrates to give 

for integration constant A. A second integration gives 

r 

Hr)=A I ^■.= AF{r,f), (2.8) 

r 

where we define a new coordinate, F, using the condition d-F := dr//. 

In principle we also must satisfy the Einstein equations (and equations for any other bulk 
fields), but instead we use the fact that we do not require more than the near-brane form for 



/(r) to side-step the effort of doing so. (See, however, [30| for many explicit solutions to the 
6D supergravity equations, including both those where the branes at r = and r = ttR have 
different properties. Many among these solutions are consistent with the near-brane forms 
being assumed here.) 

Boundary conditions and back-reaction 

We seek to eliminate the integration constants - A, f, etc. — of the bulk solution in terms 
of the physical couplings of the brane action, and this is done using the near-brane boundary 
conditions that express how the branes back-react onto the bulk Specialized to the bulk 
scalar field considered here these state 

-27rQ / - ^ = -27rQ A + gH^H + Aa^fe + /x| = , (2.9) 



where ^j, := $(0), ^'f^ := (dr^)r=o and the second equality in eq. (p.9[) uses the field equa- 
tion, eq. ( |2.7| ), as well as the form, eq. ( p.2[) , of the brane action. [One way of deriving this 
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boundary condition - for completeness, sketched in more detail in Appendix ^ - is by ex- 
cising the codimension-2 brane with a small regularizing codimension-1 cylinder (designed to 
dimensionally reduce to the above codimension-2 action when the cylinder's radius is very 
small), and using Israel junction conditions for the cylinder.] 

There are similar equations governing the near-brane form of the metric and any other 
bulk fields, but for the present purposes these just dictate how the defect angle depends on 
the value of Ub when evaluated at the classical solutions for H and Similar boundary 
conditions also apply for the brane at r = ttR, and together with eq. (|2.9|) these generically 



can be used to remove the two free integration constants in <I>(r) |26| 



The brane- localized fields must also satisfy their own classical field equations, 

§ = 0, (2.0) 
and so for x-independent H eq. ( p. 91 ) should be supplemented with^ 



H^H = ^{l_,l-g^b)- (2.11) 
This can be used to eliminate H from (^.9|), to give 

-27ra^ + A2eff$fe + /t^lcff = 0, (2.12) 
where we define the 'effective' couplings 



A2cfr :- A2 - fr-! A^leS + ^V" • (2-13) 



2A' /'i'eff-M$+ 
Divergences and classical brane renormalization 

The complication of bulk divergences enters once eq. (2.8) is used to eliminate This 



diverges logarithmically near r = due to the asymptotic limit / w r there: 

^{r) = A log(r/f) + nonsingular (as r — )• 0) . (2-14) 

Because ^ diverges logarithmically as r — )• 0, we first regularize by taking r — t- e <C i?, and 
then renormalize by allowing the brane couplings to be e-dependent in such a way that e — >• 



can be taken smoothly JT8| , 19|. Although unfamiliar in RS models, such classical divergences 
(and renormalizations) are generic to any theories with sources with two or more transverse 
dimensions (making RS models the exception, rather than the rule). Physically, these di- 
vergences arise from taking the source brane to be infinitely thin, and as such they can be 
lumped together with all of the other quantum ultraviolet (UV) effects that renormalizations 
of brane couplings would in any case have to encompass. 

With this understanding the boundary condition ( |2.12| ) becomes 

-27raA + X2es{e)AF{e, f) + (e) = , (2.15) 



*See Appendix M for the relative energetics of this solution compared with the solution H = 0. 
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and we require the singular form of the couphngs A2 cfr and cs order to determine how the 
near-brane boundary condition relates the integration constants A and r. The one condition 



that eq. ( 2.15 ) remain finite is insufficient in itself to fix the e-dependence of all couplings, but 



these are easily determined by repeating the steps of [18, 19| and demanding the finiteness 
of a few other quantities. For completeness, one way of doing this is described Appendix 
which yields the same results as earlier authors when restricted to the couplings considered 
there. 

The result for the e-dependence required of the brane couplings obtained in this way is 
simply summarized as follows, 



^-2^Pi^^^) l-2^^(e''') l-l^^C^'^-) 

X(r) = X + -i^] ; p,\lf)= , (2.16) 

2 l-^F(e,f)' ^""-^ \2^a) l-^F{e,f)' ^ ' 

where it is the renormalized ('barred') couplings that are held fixed as e — )• 0. The associated 
RG equations can be found in eqs. ( |3.66| ), below. Here f is an arbitrary renormalization 
scale, and the property F(r, r) = ensures that the bare couplings may be interpreted as the 
renormalized couplings evaluated at f = e. Given these expressions, the e-dependence of the 



effective coupling combinations appearing in eq. ( 2.15 ) are easily read off: 



Using these expressions to eliminate ^Igg and A2ofr from eq. ( |2.15| ) gives a result for A 
that is finite when e — )• 0, 

2^aA = /^^^ = -^sdH = 4^^(,) , (2.18) 

^--^ne,r) l-%SaF(r-,f) 

where the second equality gives the required relation between A and f in terms of the renor- 
malized coupling 7^1, evaluated at an arbitrary scale r. The third equality shows how this 
relation simplifies when expressed in terms of /u|,gfj defined at the renormalization scale f. 
As these expressions make clear, the dependence of the right-side of the last equality on the 
arbitrary parameter f is illusory. In what follows we will often choose the arbitrary renormal- 
ization point for which the answer is most condensed, at the expense of making a logarithmic 
dependence implicit. 

As mentioned earlier, the boundary condition at the distant brane at r = vri? imposes 
a second relation between A and r, in general fixing both and so completely fixing the bulk 
field configuration For the present purposes we leave f arbitrary, a placeholder for this 
faraway boundary condition. It is easy to track in what follows because it appears only 
through the function F, a dependence that is generically logarithmic and so quite weak. 
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With this understanding the bulk solution now becomes 



$(^) = / /fLffW J ^^^^ ^ ^2.19) 



27ra 

which, when substituted into eq. (2.11) gives the Higgs expectation as 



V 27ra 



(2.20) 



which uses eqs. ( |2.16D . This is finite as e — )• 0, as expected. The corresponding formula 
expressed in terms of renormalized couplings defined at a different scale is found simply by 
running them up or down according to ( p.l6|) . 

Eq. ( |2.20| ) shows that it is the renormalized combination fij^ that must be positive for 
H to become nonzero. Notice also that eq. ( |2.19| ) shows that it is nonzero 7l|,eff(^) that 
determines when A ^ 0, and so whether $ has a nontrivial bulk profile. Physically, this 
profile arises because /u|,gfj controls the linear couplings of $ to the brane, and having these 
nonzero precludes ^>'s near-brane derivative from vanishing. 

2.3 Higgs-Bulk mixing 

We next describe the fluctuations about this background solution, with a view towards iden- 
tifying the extent to which the H — ^ couplings cause the Higgs particle to mix with KK 
$-modes in the bulk. We only track here the mixing in the H — ^ sector, ignoring in par- 
ticular potential mixing with other bulk fields that might arise within applications where the 
scalar interacts significantly with other fields (like the metric or fiuxes) that are involved in 
extra-dimensional stabilization. This is known in particular to be a real complication when 
$ is the dilaton that arises as part of the 6D bulk supergravity multiplet [^] . 

To study fluctuations we expand the regularized action in powers of the fluctuation flelds, 

p^eff(^) A F{r,f) + (l){x,r,e)- H=^\ ° J, (2.21) 



27ra y '''''' " ^2 \v + h{x) 

where x = {x^^ and = Jl%{f) /X{f) = (246 GeV)'^ is fixed from measurements of Fermi's 
constant, Gp, in muon decay. This gives the bulk action 

Sb = - j d^x^/^ Q Q'"'d,,(l) d^c^ , (2.22) 

while the on-brane potential of the scalar sector reads 

Ub = T + Xv^h^ + Xvh^ + ^h'^ + ^ </,2(0) + gvh (j){0) + I </,(0) , (2.23) 

where all of the above are 'bare' couplings. The rest of the Standard Model looks like it 
usually does, with h acting as the usual Higgs field. 

Our goal is to compute the implications of the gv h 4>{0) term that mixes the brane and 
bulk scalar degrees of freedom, and we do so in two ways. Although somewhat redundant, 
comparing both approaches provides insight and a check on our calculations. 



- 12 - 



Dt= > + > + ...= -.0-- 

qv qv 

{hh*)= -^- + -^-'....Q..i^^- + ... = - o - 

Figure 2: Perturbative calculation of the Green's function for the h operator. Dotted lines represent 
of. and dashed lines represent ZJj^ as defined in the text. Circles are used for the gvh(l){0) vertex, 
diamonds represent a ^X2(j)^{0) vertex and Df is the A2-resunimed Green's function for the bulk field, 
Fourier transformed in the brane directions but not in the extra dimensions. 



2.3.1 Perturbative method 

A straightforward way to approach Higgs-bulk mixing [12| is to regard the terms ^X2(j)'^{0) 
and gv h(j){0) as part of the interaction lagrangian so that the unperturbed system does not 
mix brane and bulk. The implications of mixing are then found by summing all possible types 
of insertions of the mixing interactions. 

For instance, at tree level a calculation of the correlator {hh*)^ requires summing over 
the diagrams of Fig. 0. Summing these graphs gives the following momentum-space result 



{hh* 



(2.24) 



l+[i\2 +{gv)W^]Dt{0,0) 
This expression can alternatively be derived using a Schwinger-Dyson approach, as shown in 
Appendix ^ The advantage of this approach is that it does not require A2 to be perturbatively 
small. In eq. (|2.24 ), is the momentum-space propagator of h in the unperturbed theory^ 



I? 



+ 2At;2 



(2.25) 



and D'^{0, 0) is the unperturbed propagator for the bulk field, Fourier transformed only in the 
four brane directions and with both extra-dimensional positions evaluated at the Higgs-brane 
position (r = 0). 

The unperturbed bulk propagator satisfies 



-f^ + jdl + drifdr 



D^ir, 9; /, 6') = i5{r - r')5{9 - 9') 



(2.26) 



and it is useful when solving this to expand in a basis of unperturbed eigenmodes for 



ct){k,r,9) = Y,<Pni{k)Zni{r,9) , 



nl 



where 



^2Zni{r,9) 



1 

7 



Zniir,9) = -M^iZniir,9). 



(2.27) 



(2.28) 



^Notice we use e to denote the cutoff, while e governs the poles of the Feynman propagator. 
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Rotational invariance of the background further allows the separation of variables 

^nz(r,0) = P„KOe'"'/", (2.29) 

where n is an integer and the Pni are the corresponding set of radial mode functions.^ Because 
these are unperturbed modes, they do not yet 'know' about the Higgs-bulk couplings, and so 
satisfy the comparatively simple near-brane Neumann conditions corresponding to no brane 
couplings, 

{fdrPnl)r=0,.B = ^^ (2-30) 

and satisfy the traditional Sturm-Liouville bulk normalization relations 

TXR 

27ra j dr f P*iPnV = 6iv . (2.31) 



The result that follows from using this expansion in ( |2.26D is 

and, since only n = modes survive as r, r' — )■ (as is shown below), the required brane-to- 
brane propagator becomes 

Continuum limit 

These expressions become particularly simple in the large- volume limit, where \k'^E?\ S> 1. 
In this case the discrete mode spacing is very small and the mode sum is well-approximated 
by a continuum momentum integral. This continuum limit is taken explicitly in Appendix 
starting from the discrete mode sum in a simple toy model, but the near-brane result can be 
obtained more simply by solving ( |2.26| ) directly for noncompact extra dimensions. This can 
be done explicitly near the branes, where / ~ r, since eq. ( 2.26 ) becomes the equation for a 
free 2D field in cylindrical coordinates, whose solution are given in terms of Bessel functions. 
Demanding normalizability near r = for the unperturbed functions, one finds 

oo r / A \ in{9~9') 

Dt{r,e;r',e') = -z / |^ J+,2_,, ^|nH ('ZO J|nH (^O . (2.34) 

n=— oo ^ / ^ 



Since only the n = term contributes when evaluated at r = r' = 0, the near-brane limit 
becomes 

oo 



®Do not confuse these with Legendre functions, despite the notation. 
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The integral in eq. (|2.35| ) diverges logarithmically at large q, and once this is regularized with 
a cutoff A = 1/e the result becomes 



DIM 
hh 

(i.e. A;^ is timelike), in which case 



(2.36) 



The ie prescription tells us which branch of the logarithm should be used when k'^ is negative 



log (k^e^) = log {-k^e^) 



ITT . 



(2.37) 



Using this in the Higgs two-point function, and expressing the result in terms of renor- 
malized couplings gives the finite final (continuum) result 

-1 

log(A;^f^) 



{hh* 



1 



Ana 



log(A;^r^) 



(2.38) 



where f is the same arbitrary renormalization energy scale at which the renormalized couplings 
are also evaluated, and in we have used F[r, r') = log(r/r') appropriate for the large R limit. It 
is the implicit f-dependence in these couplings that cancels the explicit dependence appearing 
in the logarithms, ensuring that f does not contribute to physical quantities computed from 
{hh*)k. 

2.3.2 Direct mode diagonalization 

We next provide an alternative derivation of the Higgs two-point function, which proceeds 
more directly by calculating {hh*)k explicitly by diagonalizing the KK and Higgs modes. 
This provides a more physical interpretation for the branch cut introduced into {hh*)k by the 
Higgs-bulk couplings, in terms of Higgs mixing with KK bulk states. 

Propagation eigenstates are found by solving the field equations for h and 0, keeping track 
of the boundary conditions near the brane. In the present instance the relevant equations are 
the bulk scalar equation 



(□4 + n2)</' 

the (linearized) Higgs field equation on the brane, 

-d^d^'h + 2\v'^ h + gv 0(0) = , 
and the near-brane boundary condition for cj), 

-27rQ {fdr(t>)r=Q + 9vh + \2 0(0) = 



0. 



(2.39) 



(2.40) 



(2.41) 



As before we decompose the 6D scalar into a KK tower by expanding the solutions to 
eq. ( 2.39| ) in a basis of eigenfunctions of ^2, 



(t){k, r, 6*) = ^ •^ni>{k) Zniir, 6) , 



(2.42) 



nl 
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with 



□2Z„,(r,0) = -M2,Z„,(r,^) 



(2.43) 



We use the indices (ni) rather than (nl) here to emphasize that they run over a shghtly 
different range, with i including a value corresponding to h in the special case n = 0, in 
addition to the complete range of / for the unperturbed n = KK modes. This change only 
happens for the n = modes because only these mix nontrivially with the brane. 
As for the unperturbed case, we write 



^nK^^)=^n£(r)e^"'/", 
where n is an integer and the Vni are the radial mode functions, satisfying 



0. 



(2.44) 



(2.45) 



In terms of the 4D modes (pni and /i, the (p and h field equations, eqs. (|2.39 ) and ( p.40 ), 
are as for the unperturbed case. 



r + 2Xv' 



h + Y,[9v'Pni{0) 



(2.46) 



but Vni differs from by satisfying the near-brane boundary condition, eq. ( p.41| ), including 
the implications of the Higgs-bulk mixing. 



gvh + ^ \-2'na (fdrVne) + A2 Vni 

nl 



r=0 



(2.47) 



Using eqs. ( |2.46| ) to eliminate h and k'^ gives 
and allows ( ^.47 ) to be rewritten as a boundary condition purely for Vni'- 

{gv? 



(2.48) 



-27ra {fdrVni) + A2 



Ki - 2At;2 



0. 



(2.49) 



r=0 



What is unusual about this boundary condition is the presence of -M^^, which makes it 
mode-dependent, at least for those modes'' for which it is satisfied with Pn^(O) 7^ 0. The 
presence of the unorthodox mode-dependent near-brane boundary conditions implies the 



^The only normalizable modes for which Vni{Q) 7^ are those with n = 0, expressing conservation of 
angular momentum, as expected. 
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eigenfunctions need not be orthogonal using the usual Sturm-Liouville (or Wronskian) in- 
ner product. Instead, as shown in detail in Appendix the natural inner product adapted 
to this boundary- value problem also involves some boundary dependence. 

The resulting generalized orthonormality condition for this new inner product is 

2na JdrfV,n+ (^2_2A.2)(M,2_2A.2) = ' (2.50) 



where we collectively denote s,t = {n£}. As shown in Appendix ^, it is this modified 
boundary condition that ensures that the KK expansion of the action gives a quadratic 
action that is diagonal in the (pni- This property would be ruined in the present instance for 
an expansion using the unperturbed mode functions, (p = J2ni 't'm Pni by the term gv h (^(0) = 
gv h^^i 4>niPniiO) in the lagrangian, which causes bulk modes with n = to mix with h. In 
Appendix ^we make this discussion explicit by solving the perturbed wavefunctions subject 
to this boundary condition in an illustrative toy model. 

The complete mass eigenstates of the theory obtained including this mixing are related 
to the unperturbed states discussed earlier by a linear rotation, as follows: 

h = '^BsiPs 4>i = ^^Uis^Ps, (2.51) 

s s 

where the index i denotes {n, /} in the same manner as s denotes {n,^}. Because only n = 
modes mix with the brane, in practice tiis = Sis unless n = 0. 

The quantity Bs is of most interest, because it controls the two point function for the h 
field, which can be written 

s s * 

Eq. (I^D gives ;S, as 



_ gvVs{^) , . 

- M2 - 2Xv^ ' (2-5«i) 

and Uis is found by using 7^s(r) = Yli^is Pi{f) together with the normalization conditions, 
eqs. ( 2.31 ) and ( 2.50|) , respectively satisfied by Pi and Vg- This gives 



Y,l^*isl^^t + B:Bt = 5st, (2.54) 

i 

as well as the remaining unitarity conditions 

J^^:^, = 1, ^Z^j,S: = and ^Z^i,ZY;, = <5,,-. (2.55) 
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The continuum (large R) limit 

Although the formulae in j j2.3.2 so far make no assumption about the relative size of Xv"^ and 
the KK scale, ^ 1/^^, we pause here to display the simple result that obtains in the 

limit Xv"^ ^ fn'tiK (appropriate to large-volume models, say) for which KK mode sums are 
more usefully cast as integrals |p!o| , 

In this limit it is a good approximation to write the sum appearing in the two-point 
function as an integral 



{hh*) 



/c2 + M2 - ie 



fc2 + M2 - ie 



(2.56) 



which defines the spectral function, ph{M'^). The simplest way to obtain an expression for 
Ph is by using the explicit result for {hh*)k obtained above from the perturbative calculation. 
To read off ph{M'^) from this calculation we employ unitarity, in the form 



7:ph{M'') = Re{hh*)k 



(2.57) 



Using expression ( ^.38 ), and separating the real and imaginary parts of the logarithms, gives 
the finite result 



vr/9fe(M2 



[M2 - 2A( 



r)v^\ + vr^ 



Mf)[M2-2A(r>2] + I^ 



f=l/M 



1>2C(M2) 



where 



and 



[M2-n(M2)]2+^;4(2(J\^2) ' 

g{r)/Aa 



C(M2) 



1 



Mr) 



log(M2r2) 



+ 



Mil 

4a 



(2.58) 
(2.59) 



U{M^) := 2A(f)t;2 + 



f{r)v^\ log(M2r2^ 



Mr) 



log(M2f 



7rA2(r) 
4a 



l-(%f?)log(M2r- 



+ 



Mr) 
4a 



(2.60) 



For future reference we note here that many of the above expressions simplify if we make 
a specific choice, f = 1/M, for the arbitrary renormalization scale, so that the log(M2f2) 
terms disappear. In the present instance this leads to the simpler formulae 



and 



C(m2) 



l + (A2(f)/4a)^ 



f=l/M 



A(f) 



C(M)A2(f) 
8a 



f=l/M 



(2.61) 



(2.62) 
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Figure 3: Scattering via virtual exchange of the Higgs-bulk tower. A factor of Bs comes from the 
vertex between ips and both the initial and final state. All modes are summed over because any of the 
LPs can mediate this interaction. 



The weak-coupling limits of the spectral function are most easily seen using the second 
equality in eq. ( 2.5^ ), together with the representation 'k5{x) = \\mc^QQ/{x^ + C^). This 



gives 

PhiM"^) 6{M^ - 2Xv'^) , (2.63) 

both when g ^ and when A2 — >• oo, illustrating how the Standard Model is obtained in 
both of these limits.^ Although it is clear why this should hold when g = 0, it turns out also 
to hold when A2 — ?• 00 because the boundary conditions, ( |2.41| ), imply in this case that the 
KK modes all vanish on the brane. 

It is the two functions C and 11 that control the phenomenology of Higgs-bulk mixing, as 
we now show. 



3. Phenomenological implications 

In this section we discuss the leading sources of constraint on the Higgs-bulk mixing just 
described. We separate the effects into three types: virtual effects from exchanging KK 
modes; real effects from KK modes in the final state; and changes to the relation between the 
Higgs mass and couplings (together with associated changes to the constraints from vacuum 
stability) . 

3.1 Virtual Higgs Exchange 

Consider first a tree- level Standard Model parton process of the form X — )• /i — )• y, where an 
initial state X produces a virtual h that subsequently produces state Y, as in Fig. § The 
Standard Model amplitude for this process can be schematically written 

M'^'iX -^Y)= M{X ^ h) {hh*yk''M{h ^ Y) , (3.1) 

*We describe these limits in terms of bare couplings, since complications associated with divergences vanish 
in these limits, and the renormalized couplings A, g and A2 no longer run. 
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q (OcV) 

Figure 4: The lineshape, |(/i/i*)fcp, for the KK tower exchange with various values of (A, A2,.g,a) 
evaluated at f = (125 GeV)~^. The dotted black line is the comparison lineshape for a Standard Model 
Higgs with mass, to/j — 125 GeV and a Standard Model width of V sm — 4 MeV. (In the Standard 
Model the Higgs coupling for this mass is A = 0.1291.) The dashed blue line shows a similarly narrow 
peak (chosen to lie near 124 GeV rather than 125 GeV to avoid clutter in the figure), obtained using 
(A, A2, 5, a) — (0.127, 0, 0.0021, 0.8). The solid red line shows instead a broad peak at 125 GeV for an 
exaggerated choice of couphngs (A, A2,g,Q!) — (0.16,2,0.35,0.01). 



and so once h mixes with the KK tower this becomes 

Mix -^Y)=M{X^ h)M{h -^Y)Y, \Bsf{^s^s)k 

s 

= M{X ^h)M{h^Y){hh*)k, (3.2) 

where = = Py represents the 4-momentmxL flowing down the Higgs line. That is, 
virtual effects of Higgs-bulk mixing exchange processes are found by using (/i/i*)^ in place of 
the Standard Model Higgs propagator. 

3.1.1 Lineshape 

Implications of this modification to Higgs exchange are easiest to study in the large- volume 
limit, for which the KK sums are well-approximated by integrals. In this case we may use 
eq. ( 2.3^ ), which we rewrite in the Lorentzian form by separating the real and imaginary 
parts of the denominator 

{hh*)j^ = -i \u{-e) + e- iv'^ci-k'^)] ~^ , (3.3) 



where the functions are, not surprisingly, the same functions as in §2.3.1. 
Physical processes depend on the squared magnitude, 

\{hh*)J^ = \ , (3.4) 

which defines the resonant Higgs lineshape in the presence of mixing. A plot of this lineshape, 
and a comparison with its Standard Model counterpart, is given in Fig. ^, for several choices 
for the Higgs-bulk couplings. 
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Figure 5: The Higgs quartic coupling A required to ensure rrih — 125 GeV, as a function of the 
brane bulk mixing parameter G — gj \foL for various choices of the coupling parameter h.2 — ^^ij a (as 
listed in the legend). The flat line A2 = corresponds to the Standard Model value A = 0.1291. All 
couplings are renormalized and evaluated at a scale r = (125 GeV)~^. 




The position of the resonant maximum occurs at 

f(r) A2(f) 



n(m2) 



2A(r-) 



(4a)2 + A^(f) 



(3.5) 



r=l/mh 



where we have used the definitions of (" and 11 in eqs. (l2^ ) and (|2!6^ ). This suggests the 
definition 

m2:=2Aefr(rV with 2K^{f):=2\(f)-^^^^^ (3.6) 

f=\lm,K (4a)^ + A2(r) 

It is the value of the renormalized couplings at the scale f = l/nih that is relevant to 
many physical quantities, such as the condition AefT(f = l/rrih) = 0.1291 that is required to 
ensure a Higgs mass of 125 GeV. Because of this, in what follows a barred coupling without 
a specified renormalization scale is understood to be evaluated at l/nih- 



A = A(f) 



g = g{r) 



and so on . 



(3.7) 



Notice also that the condition m/j = 125 GeV imposes only a single relation amongst the 
couplings A, A2 and g, rather than fixing A completely, as it would have done in the Standard 
Model. Fig. || plots the value predicted for A as a function of g for various choices A2 , and we 
assume in what follows that A is fixed in this way. 

We see that the lineshape can resemble a single Higgs resonance despite its containing a 
sum over many KK states. Its width at its maximum is 

g'^v'^/4a 



1+ (A2/4a)^ 



(3.8) 
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This width is related by unitarity to the rate for invisible processes where the KK modes 
escape invisibly into the bulk, carrying with them missing energy. When Fg is sufficiently 
small, such as in the g limit, its role in the unitarity argument is eventually replaced 
by the Standard Model Higgs decay width, Tsm, corresponding to the imaginary part of 
the usual Standard Model Higgs vacuum-polarization graphs. This suggests the definition 
Th = Fg + Ts]^f. However, because the Standard Model contribution is a loop effect it should 
only be kept in the special case where it dominates at the peak of a narrow resonance. 
Given these considerations, we write the corrected Green's function 

{hh% = -i [n(-fe2) + fe2 _ iy^(^(-k^) _ imhTsM] , (3.9) 

so that ^ 

" m-k^) + P]2 + ruh [vKi-k^) + rrihTsM] ' ^^'^^^ 

3.1.2 Invisible width 

In essence, mixing with the bulk introduces a new invisible channel into Higgs reactions 
while leaving unchanged the relative strength of all of the visible h couplings to other Stan- 
dard Model particles. The exchange of the KK tower (instead of just the Higgs) suppresses 
the overall rates for observable Higgs-mediated processes, while preserving their relative fre- 
quency. The success of the Standard Model description of the resonance at 125 GeV, as seen 
by both CMS and ATLAS, therefore provides an immediate constraint on such mixing, in 
much the same way as it constrains a more conventional Higgs invisible width. 

In the present context the nature of this constraint is easiest to see within the narrow- 
width approximation, which applies in the phenomenologically most relevant case where 
m/j 3> -I- FsM- In this case the resonant h autocorrelation becomes 

rrih [i s + J- sm\ 

which neglects a factor [l — H'(m^)] ^ . This factor can be dropped because 

and so is similar in size to other contributions that have been neglected. 

For comparison, the Standard Model Higgs distribution in the same narrow-width limit 
reads 

|(M*),|L « Sik'+ml) , (3.13) 

'l''h^ SM 

with mh being the physical mass of the Higgs. In this limit the momentum dependence of the 
two results is identical, and the exchange of the KK tower just provides an overall suppression 
to the rates for Higgs-mediated processes by the factor 

n = , (3.14) 
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Figure 6: Bounds on A2 = A2/Q! and G = g/-/a from constraints on the invisible branching ratio of 
a Standard Model Higgs. The shaded region is excluded, for two sets of assumptions (described in the 
main text). Couplings are evaluated at a scale f = (125 GeV)~^. 



relative to the Standard Model. 

The bound that follows for IZ can be inferred using the results of extant global fits to 
the LHC data that were performed to constrain the branching ratio into invisible decays 
of an otherwise Standard Model-like Higgs. The corresponding narrow-width suppression 
for a conventional invisible decay width would be TZ\av = '^sm/(^sm + ^im) = 1 — Bmv, 
where i?inv is the branching fraction into invisible decays. Recent fits give an upper bound 



B\nv ^ (0.3 — 0.64) at 95% CL |J, 35], where the range depends on precisely the priors 

used when performing the fit. (Ref. [^3| finds Bi^v < 0.64 using 15 signals from the Tevatron, 
Atlas and CMS, while ref. p2| finds the slightly stronger limit of i?inv < 0.4 at 95% CL using 



a different suite of 16 observables. The bound from ref. |35] is similar. Ref. [34| finds the 
strongest bound, i?inv < 0.2 using a smaller set of Higgs signals that are argued to be more 
sensitive.) 

Taking the most conservative of these limits, we infer the constraint 

7^>0.36. (3.15) 

The corresponding constraint in the g — X2 plane is plotted in Fig. ^, which also shows the 
result obtained from the more aggressive constraints. This plot shows that these global Higgs 
fits imply a conservative limit 

^ < 0.007 for A2 = ; (3.16) 

a 



with weaker constraints on g as A2 increases. (This weakening of the g constraint with large 
A2 is a general consequence of the decoupling of brane and bulk in the A2 — )• 00 limit.) 
For comparison, Ref. [|l^] considered the phenomenological implications of the cubic vertex 
^gh^(j){0) but neglected Higgs-bulk mixing and A2. They found that g = 0.18 was accessible 
at the LHC with 100 fb"^ at 14 TeV in the /i 7 7 (/> channel. 
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Figure 7: The Feynman graph corresponding to the one-loop correction to the muon anomalous 
magnetic moment in the Higgs-bulk mixing scenario. 



Notice that Tg depends on the two variables g and A2 only through the combination 
C{mf^), and this is generally true (once ruh is fixed) for any observables for which the narrow- 
width approximation is justified. Whenever this is true it is more useful to quote the constraint 
directly on C{''^h)i giving 

C(m2) < 5 X 10-5 . (3.17) 



3.1.3 Low-energy bounds 

The effects of other virtual contributions of KK modes can be similarly computed given 
an expression for (hh*). Important among these are constraints from low-energy precision 
measurements, such as the anomalous magnetic moment of the muon. By way of example we 
study this here for the large- volume limit where (hh*) has a simple closed-form expression. 

This section uses this calculation to conclude that these bounds are negligibly weak, in 
agreement with standard intuition for the large-volume case. This standard intuition starts 
from the observation that each KK mode couples with gravitational strength, and so it is 
only the enormous phase space for KK modes that can compensate for this suppression 
[|lC|]. Bounds on extra dimensions from low-energy observables are usually weaker than those 
from colliders because, for low-energy processes, only sub-TeV modes contribute, making the 
phase-space compensation incomplete and so insufficient to produce an observable result. 
Astrophysical energy-loss bounds on extra dimensions are an important exception to this 



intuition |3^, |3^, and we return to these below. 

Fig. ^ displays the Feynman graphs whose evaluation gives the KK mode contribution 
to a fermion's anomalous magnetic moment. We estimate that the difference between this 
graph and the corresponding graph for a Standard Model Higgs is quite small, as it should 
be if we consider perturbing in small g, and brane-bulk mixing effects are negligible. 

To see this, consider first the Higgs contribution from the analogous graph in the Standard 
Model. This evaluates to an anomalous fermion magnetic moment of size 

Sn'^mf^ J (1 — x) -I- x^(m^/m?) ' 
■' 



where is the fermion's mass and yf := rrif/v is its Higgs Yukawa couphng. The regime of 
practical interest is mj ^ mh and, because the integral diverges logarithmically near x = 1 
when mj = 0, the answer in this limit is dominated by 




(3.19) 



where the numerical values assume the fermion is the muon. This is negligible in comparison 
to both the electroweak boson contributions and the experimental precision, which are of 



order a few 100 x 10"^^ |0|. 



Repeating this exercise for the graph in Fig. ^ using the large-volume expression for 
{hh*)j^ allows us to estimate the difference between the Standard Model Higgs contribution 
ah and the contribution from the whole KK tower a^. We do not use {hh*)'f^ with the inclusion 
of the Standard Model Higgs width because it is a higher-loop effect, and for simplicity we 
assume A2 = 0. Therefore, accounting for the new propagator in the graph is accomplished 
by making the replacement 

1 1 1 (fv'^\ log(A:2f2) 



+ 



A:2 + 2A7;2 ^2 + 2A(f)t;2 + (^gj) log(/fc2f2) k'^ + 2\{f)v'^ V 4^" 7 [fc2 + 2A(f)t.2] ^ ' 

(3.20) 

where we write g{r) = g because it does not run in the limit of vanishing A2. 

Although the above propagator is only one term in the graph, which also contains a loop 
integral over k"^ , in order of magnitude, we estimate that the relative difference between the 
anomalous moment in the Higgs-bulk scenario and in the Standard Model is 

- (3.21) 



ah 



SvraA 



and so can be ignored given the strong constraints already found for g. 
3.2 Invisible final states 

This section considers the implications of Higgs-bulk mixing for observables with Higgs-bulk 
final states, rather than simply as intermediate states. We first argue that these states 
are invisible and then relate their production rates in various channels to the analogous 
Standard Model Higgs production rate. Invisible states appear as missing energy at high- 
energy colliders, and we discuss Higgs-bulk mixing signals at both LEP and the LHC. We 
then turn to low-energy constraints on missing energy in astrophysical systems. We conclude 
that LEP provides weak bounds on the Higgs-bulk couplings, but astrophysical bounds are 
considerably stronger. We also provide a preliminary estimate of the reach at the LHC, beyond 
the constraints already discussed coming from the suppression of the rates for producing 
Standard Model particles. 
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3.2.1 Invisible-state production rates 

Consider a process whereby the tower of states is produced through their overlap with 
h. For weak couplings the physics of this process resembles the physics of sterile neutrinos, 
in that production and detection of the produced state occurs only because of the mixing 
with a weakly-interacting Standard Model particle. For weak couplings invisible processes 
are therefore described by requiring the amplitude for producing a state orthogonal to h; that 
is, for being in one of the 'flavour' eigenstates, (pi rather than a propagation eigenstate, ips- 
We therefore consider the squared amplitude to produce a (pi final state 

\M{X ^ y <A.)|' « \Ms,,{X ^Yh)\l^ \M{h ^ (Pi)\l^ , (3.22) 

where the first factor is the Standard Model result for Higgs production, but evaluated using 
P(f), the momentum of the final state (pi, which in practice amounts to replacing m/i — >• Mj. 
The second factor can be understood in two ways. Formally, M{h — t- (pi) is given in terms 
of {h(p*)k by amputating the final unperturbed (pi propagator and putting the correlation on 
shell (as usual) to obtain an amplitude with (pi in the final state: 

M{h ^ (P,) = Amp [{h(P*{0))k=pJ := {h(P*mk=p, [Dt^y . (3.23) 
As shown in Appendix Q the amputated {h(p*)k correlation is given by 

ig[f)v 



Amp [{h(P*{Q)), 



l + i(A2(r)/4a) 



{hh*), (3.24) 



where the renormalization point r = 1/Mi is chosen for notational convenience. Alternatively, 
the factorization in eq. (|]2|) and the expression in eq. ( p. 24]) can be derived by directly 
summing the graphs of Fig. ^ However, the virtue of using a Schwinger-Dyson equation over 
the graphical methods is that it emphasizes that the result need not assume A2 is small. 



Using this in eq. (3.22) then gives the general result 



\M{X ^ y 0.)|2 « Aav\{Mf)\{hh*)k=pf \Msm{X ^ Y h)t , (3.25) 



where {h h*)k=p^ is given explicitly in terms of the unperturbed h and (p propagators in earlier 
sections (c./. eq. (|2.24| )). 
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In practical applications, it is the sum over modes' squared amplitudes (possibly weighting 
some other function O) that is relevant to a physical observable. That is, 

O = Y,\MX ^Y(j)i)\'^0{M^). (3.26) 

i 

This sum over bulk states, and the expression for {hh*)k=p^, both simplify considerably in 
the large- volume limit. In this limit we can also account for the on-resonance effects of decays 
into Standard Model particles, by using the large- volume expression for {hh*)'j^^p^ in place of 
{hh*) i^^p^. Then, the observable can be written in the simple resonant form 

O « y dA/2 T(Af2) 0(m2) \Msm{X ^Yh)\l^ , (3.27) 

where 

^ [M2 - n(M2)]2 + [v2({M^) + TUhTsM? ■ ^ ' 

In taking this continuum limit we use the fact that each mode has only a gravitational- 

— 1/2 

strength coupling to the brane, due to the proportionality of each normalized mode to V2 , 
where V2 is the extra-dimensional volume. Consequently, it is only the enormous phase space 
available at high energies that can compensate for the extremely feeble coupling of each mode, 
implying that it is only the density of states of the high-energy modes that is important. But 
the density of these modes does not depend on the details of the shape of the extra dimensions 
(unlike the density of states for the lowest-lying modes), which allows us to use the flat-space 
result appropriate to an extra-dimensional torus, 

' ^I{m!) « / /(p2) = / ^ I{M') , (3.29) 



V2 ^ J (27ra)2 ' J Ana 

i 

even for applications to more complicated geometries like spheres and rugby balls. 
3.2.2 Missing energy at LEP 

If the Higgs mixes significantly with invisible light states then these should have been produced 
at LEP, leading to a constraint on Higgs-bulk mixing. A convenient way to obtain this 
constraint is to use a particular search performed by the ALEPH, DELPHI, L3 and OPAL 
experiments at LEP II. These experiments have sought evidence for Z boson production in 
association with a Higgs that decays 100% invisibly while being produced with a Standard 
Model rate |41|. A combined analysis of each experiment's ^/s = (200 — 209) GeV data has 
been used to place a lower bound on the mass of such a Higgs as mh ^ 114.4 GeV at 95% 
CL. 

This is a convenient search for the present purposes for two reasons. First, the resulting 
bounds usefully constrain the cross section at these energies for generic electron positron 
annihilation into a Z boson plus missing energy, since it gives the same signal; that is 

fTexp(e+ e- ^Z$^)< asM{e+ e- ^ Z h) , (3.30) 

771^=114 GeV 
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where both sides are evaluated at ^/s = 209 GeV. This is true so long as the selection 
efficiencies are similar for the new process and the Higgs process. Second, as we expect from 



the discussion in § 3.2.1 , the invisible cross section predicted from Higgs-bulk mixing shares 
enough of the features of the Standard Model cross section to allow a simple inference of the 
constraints. 

Any of the (pi states produced at LEP would appear as missing energy, so that the total 
missing energy cross section is the sum of the individual cross sections 



(3.31) 



The total missing energy cross section is therefore a weighted sum of squared (pi production 



amplitudes. Using the results of §3.2.1, we can write the cross section for the missing energy 
process in the form of an integral 



M2 

-"-'max 



a(e^ 



dM- 



;T(M2 



(l-KflE, (27r)32S<^ 



(3.32) 

where the upper bound on integration Mmax = \/s — ruz with ruz as the mass of the Z bo- 
son. This reflects the fact that only these modes are kinematically accessible. Additionally, 
|A^5mIp^ = \M.sM{e'^e~ — ^ is the Standard Model Higgsstrahlung amplitude appropri- 

ately spin-summed/averaged, with the subscript, p^, reminding us that it is evaluated using 
the (pi final-state four-momentum p^ = {E^,p^) where E'^ = + M^. Here and in the 
following, is the usual initial-state-dependent flux factor associated with a cross section, 
Px is the four- momentum of this intial state and pz = {Ez, Pz) is the four- momentum of the 
Z boson. 

The Standard Model Higgsstrahlung cross section can be written 



o"sA/(e"^e Z/i) = ^ J 



d3 



Pz 



{2T:f2Ez {2-Kf2Eh 



\MsM\lA27i)H^{px-Pz+Ph) 



(3.33) 



with Ef^ = pf^ + mj: where ph = {E^, ph), and sm\'^ is evaluated using the Higgs momentum 
Ph- We recognize the same expression in the missing energy cross section with the replacement 
nifi — )• M. This allows the missing energy cross section to be rewritten as follows 



A/2 



(Tee 



Z$r)= j dM^T{M^)asMie^ e- ^ Zh) 



mh=M 



(3.34) 



This expresses our missing-energy prediction in terms of the well-known Standard Model 
Higgsstrahlung cross section asMi whose m/j- and y^-dependence is given by a standard 



result ^ 



(Ts„(e+ e- ^Zh)(x ^J{ml) [j(m^) + 12M|/s] , 



(3.35) 
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A2 



Figure 9: Constraints from LEP expressed in the A2-G plane (with G = g/y/a and A — \2/a) 
from LEP constraints. The small shaded region is excluded. All couplings are evaluated at a scale 
f = (125 GeV)-i 



where J{m?) := (1 — m^/s — Mljsf' 
Using this, we may write the bound 



^m^M'^js^ is the two-body phase-space function. 



cr(e+e < crexp(e"^e Z$j,) < asMie~^e ZH) 

in the form 



mh = 114GeV 



(3.36) 



dM^ T(m2)VJ(M2) [J(m2) + 12Af|/s] < ^J(m2) [J{ml) + 12Af|/s] 



mh = 114GeV 



(3.37) 

where all other pre-factors in the cross-section formula cancel. The integral can be evaluated 
numerically to determine the allowed region of parameter space, with the result plotted in 
the 5 — A2 plane in Fig. The result excludes an island in parameter space, and at first 
sight it might appear surprising that the bounds get weaker at large g as well as at small g. 
The weakening for large g is a consequence of the resonant shape of T{M'^) together with 
the proximity of \/s — niz = 118 GeV to the resonance's maximum, because T oc \/{C,v'^) 
when (^f 2 is bigger than both — Il{M'^) and uihT sm- (A similar thing happens for bounds 
obtained for other observables dominated by the resonance, though the weakening of the 
constraints at large g for these occurs for g too large to justify our approximations, and so is 
not shown in the plots.) 

The bounds obtained are clearly weaker than the LHC bounds on the Higgs invisible 
width discussed above (and the astrophysical bounds discussed below). So, although the issue 
of selection efficiency was ignored in arriving at these constraints, their weakness illustrates 
that a very large discepancy in the two signals' selection efficiencies, that furthermore favours 
Higgs-bulk events, would have to exist for LEP bounds to become significant. 
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3.2.3 Missing energy at the LHC 

The dominant production mechanisms for a 125 GeV Higgs at the LHC are gluon fusion 
gg ^ h through a top-quark loop and vector-boson fusion qq — )• qqh. Additionally, the 
Higgsstrahlung process qq — )• Vh provides a clean signal at the cost of a reduced cross section 
with respect to the other production channels. 

These processes give rise to three different missing-energy signals for an invisible Higgs. 
At non-leading order, additional jets can radiate from coloured particle lines in gluon fusion, 
resulting in j and jj final states, the second of which is also attainable at leading order 
through vector boson fusion and Higgsstrahlung with a Z decaying hadronically. Alterna- 
tively, a Higgsstrahlung Z can decay leptonically, giving a signal i'^i" 

Missing energy rates from Higgs-bulk mixing at the LHC are also determined from the 
rate to produce 4'i final states. The arguments in up to eq. ( |334D can be repeated 

for any of Higgs production processes at LHC, at least at parton-level, so that the missing 
energy, parton level cross section, cr, is 



M<\/I—mY 



a{X = j dM^T{M'^)asM{X ^Yh) 



=M 



(3.38) 



where s is the parton-level center of mass energy squared, and my is a placeholder for the 
region of M phase space denied to (pi by the presence of the other final state particles. For 
example, if the final state particles associated with (pi are all massless, = 0, and for 
Higsstrahlung my = ruz as in § |3.2.2 . 

An important difference between the LHC rate and the rate at LEP is that the LHC is 
sufficiently energetic to probe the peak of T(M^) that lies at rrih- This allows us to employ 
the narrow width approximation for rrih ^ Tg + r^jv/) in which 



ml) , 



(3.39) 



where we have used the fact that v'^C,{mh) = rrihTg. This simplifies the parton-level missing 
energy cross section 



Ts + r 



SM 



asMiX^Yh). 



(3.40) 



One might worry that there are collisions where the initial state partons have small 
momentum fractions xi,2, resulting in a small center of mass energy \/I = ^xiX2S. If is 
too small, then the integral over will not saturate the delta function 5(M^— m|). However, 
the resulting step function Q{'mh+mY — also appears in the Standard Model cross section 
where it encodes the threshold above which the Standard Model Higgs can be produced, and 
so we refrain from rewriting it. This calculation reveals that, in the narrow width limit, we 
expect the parton-level rate of producing an invisible Higgs to be suppressed relative to the 
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Standard Model cross section by a factor of the invisible branching ratio B = T b/ {^b + '^sm), 
as in other models with additional invisible Higgs decay channels. 

The trivial mass and energy dependence of this factor allow us to convolve the parton 
distribution functions without complications. Furthermore, since it is only the modes near 
M ~ rail that contribute to the cross section in the narrow-width approximation, kinematic 
cuts will apply equally to single Higgs cross section and the Higgs-bulk cross section. There- 
fore, for a given production mechanism, the total missing energy proton-proton cross section 
at the LHC can be written 

a{pp ^Y$^) = Bx asuipp Y h) , (3.41) 

where B = Tg/{Tg + Tsm) as before, and UsMipp -^Y h) \s the Standard Model production 
cross section at the LHC. 

New-physics constraints and discovery estimates for a Higgs with invisible decay channels, 
specific to a given production channel are cast in terms of the quantity 



?^:= i?inv, (3.42) 

where Gbsm is the production cross section of the proposed invisible Higgs, Usm is the Standard 
Model Higgs production cross section, and i?inv is the invisible branching ratio predicted by 
new physics. Although these constraints envision a single Higgs state with new decay (and 
possibly production) modes, they apply equally well to Higgs-bulk production in the narrow 
width approximation, where we have assM = ^^sm and Bjnv = B and so = B. The 
constraints and reach of LHC invisible-Higgs searches are readily translated into constraints 
on Higgs-bulk mixing (as in § 3.1.2| .) 



Using the ATLAS monojet search sX ^/s = 1 TeV with 1 fb ^ Q, the authors of Ref. @ 
were able to bound ^ 20 for a 125 GeV Higgs-like particle. The updated CMS monjet 



data [46 1 with 4.7 fb -"^ at 7 TeV was subsequently used |47| to tighten this constraint to 
< 1.3 at 95% C.L. and it was argued that the bound could be increased to < 0.9 with 15 
fb^i at 8 TeV. It is only once < 1 becomes possible that measurements become sensitive 
to Higgs-bulk mixing. However, part of this region is already accessible (and ruled out) by 
constraints from the LHC global fits discussed earlier, since they require B\rsv ^ 0.64. 

Nonetheless, higher integrated luminosity, and different channels, will allow Higgs-bulk 



mixing to be further probed. Ref. |48] estimates that the LHC with 20 fb~^ at 7 and 8 TeV 
can exclude invisible rates down to ~ 0.4 at 95% C.L. via the dijets plus missing energy 
signal. Finally, it was estimated that ~ 0.25 would be probed at 5a by the LHC with 300 
fb^^ at 14 TeV [^], also via dijets plus missing energy. Both of these studies also considered 
the LHC sensitivity to Higgsstrahlung process, where the Higgs is invisible and the Z decays 
leptonically, but found it to be a weaker probe of Monojet searches were also concluded 



to be a weaker probe than vector boson fusion in ref. [48|. 
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3.2.4 Astrophysical constraints 

The new, invisible (/ij states can carry energy away from stars and supernovae. If the new 
emission process were too efficient, then it would confiict with the current understanding of 



stellar evolution. On this basis, it has been argued very generally |36] that any new energy 
loss channels in the sun must not exceed the solar energy loss rate Ssun ~ 2 ergs^^g^^ and 
new channels in red giants and horizontal-branch stars cannot exceed £^0 ~ 10 ergs~^g~^. 
Additionally, neutrino observations from SN1987a suggest that a new channel must not release 
energy at a rate exceeding the neutrino rate £sn ~ 10^^ ergs~^g~^ during core collapse. 

In this section, we consider a variety of (j) emission processes in these settings and deter- 
mine the strength of the corresponding bounds on Higgs-bulk mixing. Since the goal of this 
section is to estimate constraints, we make various simplifying assumptions and neglect the 



effects of dense media, interference effects from multiple scattering |36, 50 1 and the possibili- 
ties of trapping the KK modes inside of an astrophysical object, or of their decaying before 
exiting the astrophysical medium. 

Electron-positron annihilation in supernovae 

The temperature in the core of SN1987a Tg^v = (20 — 60) MeV was high enough to produce 
electron-positron pairs. These pairs could have subsequently annihilated into an h state that 
then mixed over to a (pi mode that carried energy away from the core. We calculate the 
associated energy loss rate for this simple process as an example to guide our discussion of 
astrophysics constraints on Higgs-bulk mixing. 

The energy loss rate for this process is given by the sum of the individual (pi emission 
rates 

£ = Y,^i- (3-43) 

i 

where the emission rate to a single state (pi is given by 

In the above equation p = 3 x lO^^g cm~^ ^ ^qcd is the density of the supernova core, 
|Alp = |A^(e^e~ — )■ (pi)]'^ is the spin-summed amplitude for the annihilation process, P(i,2) 
is the momentum of the initial state electron (positron) with energy 2) ^-iid p^ is the 
momentum of the outgoing KK mode with energy E^p. The occupation numbers /(i,2) 
given by the Fermi-Dirac distribution for a relativistic electron (positron) 

/(i,2)= J.,^ , (3.45) 

where /ig ~ 345 MeV is the chemical potential for the electron in the supernova core. Because 
the final state KK mode is assumed to escape, there is also no need for a Bose-Einstein final- 
state factor for it in the energy-loss rate. Eqs. ( 3.43| ) and ( [3.44 ) show that the total energy 
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loss rate is a weighted sum of squared (pi production amplitudes, so the results of § 3.2.1 allow 
us to write the total energy loss rate in the form of an integral over 



oo 





(3.46) 

where |A^smP = |-^sA/(e^e^ — t- is the spin-summed Standard Model Higgs amplitude 
squared. 

As in § 3.2.2 and § 3.2.3| , we could proceed by noting that the energy-loss rate can be 
written 

oo 

£ = [ dM^T(M^)£sMie^e- ^ h) , (3.47) 

J mh=M 


where £sm is the Standard Model Higgs emission rate, and it should be evaluated for a Higgs 
with mass m/j = M. However, since £sm is unknown, it is more straightforward to determine 
the total energy- loss rate via eq. ( |3.46| ), a task to which we now turn. 
The spin-summed, squared, Standard Model amplitude is 

\MsM(e+e~ ^ /i)|2 = 4{y^fE^E2{l - cos^ia) , (3.48) 

where is the Standard Model Higgs-electron Yukawa coupling, 9i2 is the angle between 
the momentum vectors of the electron and positron, and the mass of the electron is neglected 
because Tsn ^ rne- Using standard techniques to integrate phase space, we find 

oo 

£ = A^^- / d^i dE2 dM^ T(m2) /i/2(Si + E2) 9(4^1^2 - M^) , (3.49) 
oAtt-^P J 


where is a step function encoding the threshold above which a (pi mode with mass M can 
be produced. This integral must be integrated numerically using the expression for T(M^) 



in eq. (|J|). 

The low-energy form of T(M^) simplifies greatly if we specialize to A2 = 0, and consider 
the limit justified by existing constraints. First, eq. ( |2.59| ) gives (^{M"^) = g'^/Aa as 

constant. Additionally, assuming g 1 allow us to approximate eq. ( 2.60]) as n(M^) = m|. 
Finally, since T^j^ <C tti^ we find that in these three limits X(A^^) is a constant 

T(M2)«^!!^, (3.50) 

where this approximation is better than 1% for M < 1 GeV and g/a < 0.007. This simplifies 
the total energy-loss rate to 

00 
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where x = Ei/T and y = E2/T. This gives 



S « {^-^ (4.4 X lO^^ erg s-^g-^) T^oA^ijT) , (3.52) 

where T20 = T/(20 MeV) and I is the integral factor in ( |3.51| ). It is exponentially sensitive to 
fj,e/T and ranges from 2 x 10~^ to 4.4 as temperature ranges from 20 to 60 MeV. However, even 
the most severe bound on Higgs-bulk mixing from electron-positron annihilation in SN1987a, 
which is found by assuming Tgjv = 60 MeV, only constrains gj \/a <, 0.05, which is already 
ruled out (for A2 = 0) by the LHC global-fit constraints in Fig. 

Stellar processes 

Rather than calculating additional subdominant constraints on Higgs-bulk mixing, we next 
use existing constraints on the Yukawa coupling of a single light scalar ■0 to estimate whether 
a given emission process will result in a significant bound. For example, the strongest stellar 
bound on the Yukawa coupling of a light scalar to nucleons, y% < 4.3 x 10~^^, comes from 
the Compton process A 7 — )• Atp in red giants, where A = {^}ie,p} and the temperature 



of red giants is T^q ~ 10 keV 51]. The bound is derived by assuming that the light 
scalar couples only to to nucleons, which is a good approximation to the Higgs-bulk mixing 
scenario since the (pi states couple to the stellar constituents through the h state, and the 
Standard Model Higgs-nucleon Yukawa coupling is much larger than the electron Yukawa 
couphng ^ (340 MeV/v) > ||]. 

To estimate the energy loss rate for the Compton process ^ A(pm Higgs-bulk mixing, 



we write it the same form as eq. ( 3.47 ) 



00 

£ = [ dM^T{M'^)£sMiA-f ^ Ah) , (3.53) 



where £sm is the Standard Model Higgs emission rate. In order of magnitude, we assume the 
Standard Model Higgs emission rate is related to the emission rate of a single light scalar ip 
as follows ^ 

£s^fiAJ^Ah) £^{A^^Ai;)e{T^a-M), (3.54) 

where £^ is the emission rate of a light scalar. We include a step function, 0, to account 
for the fact that only Higgses of mass M < T will have appreciable production rates. This 
follows from the fact that A is much heavier than the energy of the photon, so we can neglect 
its recoil and assume that the outgoing (pi mode has the energy of the incoming photon 

~ E^. The emission rate for heavier (pi modes is therefore suppressed by a Boltzmann 
factor e"-^^/^ ^ g"*^/^ <C 1. 

Using this approximation, and taking the A2 = 0, small g, low-energy limits for T(M^) 



in eq. ( 3.50 ) then gives an order of magnitude estimate for the Higgs-bulk energy loss rate in 
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terms of the ip emission rate 



This allows the bound on to be translated into a bound on Higgs-bulk mixing (with A2 = 0) 



< (4.3 X 10-") , (3.56) 



from which it is estimated that even the tightest constraint g/y/a < 0.8 from stellar physics 
is subdominant to constraints from LHC global fits in Fig For completeness, we note that 
there are similar constraints on the Yukawa coupling of a light scalar to nucleons from the 
Compton process in the sun, and the bremsstrahlung process Ae ^ Aeip in the sun and red 
giants [^], but none of these are estimated to give improved bounds on Higgs-bulk mixing. 

Photon annihilation 

Photon annihilation is a relevant process in both both stars and supernovae. However, the 
Standard Model Higgs couples to photons through W and heavy quark loops, and so must the 
KK modes. For processes like photon annihilation 77 — )• h, in which the photons and Higgs 
are all on shell, the effect of these loops is well captured by the following effective Lagrangian 

C = -yc^F.^F'"'h with c^^ — . (3.57) 



We estimate that 



a2 r^' 



\^sMh7^h)\'^^-^^j IMsuie-^ e- ^ hr , (3.58) 

from which it follows that energy loss rate from photon annihilation is less than the rate 
from electron-positron annihilation, in stars and supernovae, and therefore gives a negligible 
bound. 

Nucleon-bulk bremsstrahlung in supernovae 

In addition to the bounds from the sun and red giants, there is a similar bound on the 
nuleon Yukawa coupling of a light scalar from SN1987a data < 4 x 10^^^ that follows from 
considering the nucleon bremsstrahlung process N N — )• N N [^3| and assuming Tsn = 60 
MeV. This can be used to estimate the constraint coming from the energy loss process N N ^ 
N N (j). Although this upper limit is similar to the stellar bounds on the same coupling, it 
provides a much more stringent constraint on Higgs-bulk mixing because the supernova core 
is much hotter than red giants, thereby enhancing the Higgs-bulk energy loss rate by a factor 
'^sn/'^rg ~ 10^- Generalizing eq. ( 3.55| ) to bulk- nucleon bremsstrahlung, we expected a 



constraint gj^fa < 10 ^ if the supernova temperature is assumed to be Tg^r = 60 MeV, 
which is dominant to all other constraints discussed so far. 
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Figure 10: The four uncrossed graphs for Higgs-nuclcon bremsstrahlung. 



In order to confirm this estimate, we calculate the energy loss rate in the one pion 
exchange approximation, using the following interactions 

C = -ig^MNNl5N - y'^.hNN , (3.59) 

where N represents the neutron, vr" the neutral pion, and ^Trivjv ~ 13 phenomenologically. This 
approximation is known to overestimate the rate of axion and KK graviton bremsstrahlung, 
where it can be tested against model-independent calculations [54|. Since model- independent 



methods are not applicable to scalar radiation via a Yukawa coupling [55|, the emission 
rate and associated bound are just crude estimates. We are content to further simplify the 
calculation by neglecting bremsstrahlung from protons, because the proton fraction in the 
supernova core is small, and by neglecting the Higgs-pion coupling, which is a dimension-5 
operator. We also neglect a suppression due to multiple scattering effects [|^] that is expected 
to reduce the emission rates of scalars in nucleon bremsstrahlung by as much as a factor of 5 

The total energy loss rate for this process is given by 

oo 



j=l,..A 



x\Msm\;^ {2TTf5\pi +P2-P3-p^-k), (3.60) 

where |A^sjv/P = \-Msm{N N — )• N N h)\'^ is the spin-summed Standard Model Higgs ampli- 
tude squared, pj represent the neutron four-momenta, and everything else maintains its old 
definition, except that it is now appropriate to use Maxwell-Boltzman distributions for the 
non-relativistic neutrons with number density 

f, = !llL(J^\"\,,J._d_], (3.61) 

and in the phase space measures Ej — )• m^v, where mjv is the nucleon mass. 

We proceed by calculating the spin-summed Standard Model Higgs amplitude |A^sa/P- 
There are 8 relevant diagrams for this process, four of which can be found in Fig. OQ. The 
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other four correspond to the crossed analogues of those hsted. In the nonrelativistic hmit 
IPjI ^ "rripf, and we assume the relevant KK modes are much lighter than the neutron M <^ 
rriff. In this limit, the squared, summed matrix element reads 



\M 



SM 



2 _ sjiAfsi, _ 



, (3.62) 



where k = p2 — P4, 1 = P2 — P3- The individual terms in the square brackets are C(l) since 
|k|, |1| ~ \/'m,^T ~ TTT-TT so that we approximate the matrix element as a constant ||36t] 



l-MsMp - ^(^^^y^^ . (3.63) 



"AT 



We evaluate the phase space integral in the non-degenerate limit. Although the neutrons 
in a supernova core are somewhat degenerate, in more detailed calculations of scalar-nucleon 
bremsstrahlung [55j, the non-degenerate limit approximated the full emission rate well for 



moderate supernova temperatures T^n ^ 20 MeV. In this limit we neglect the blocking factors 
1 — /s 4 ~ 1 ) and after integrating over phase space we find 



^ ^ < {y^)^gt^^T-' ^ I ^^^^2) ^g-y-x / 2 ^ ^ 1/2 / 2 _ ^ 1/2 ^ 

1287r7/2M„/ p ^ 

_ nliVn? Q^nnT'"'^ 1, f ,^^r.n^2.^^-xn , ^„/.^l/2/^2 n1/2 



(x^ - z) 



128 7rV2My^ 



P 



/d./d.T(.^V.-(l+W4)"^.^-.)"^ (3.04) 



where x = E^p/T, z = M"^ /T"^ and the approximation made in the second line is good to within 



2.2% [36|. This integral can be evaluated numerically, and the allowed region of parameter 



space from SN1987a constraints is plotted in Fig. 11. 

This somewhat crude calculation and the associated bounds verify that bulk-nucleon 
bremsstrahlung in SN1987a is likely the strongest constraint on Higgs-bulk mixing and, for 
Tgjv = 60 MeV, gives gj^fa < when A2 = 0, as estimated. However, since we ne- 

glected many effects that might suppress energy loss rate, and SN1987a bounds from nucleon 
bremsstrahlung are at best order of magnitude estimates, we assume the Tsn = 20 MeV 
bounds in order to be conservative. 



3.3 Mass-coupling relations and vacuum stability 

A great virtue of the Standard Model Higgs is the tight connection between the strength of 
its coupling to a particle and that particle's mass. In particular, knowledge of the mass of 
the Higgs itself reveals the strength of its self-coupling. Although it is not yet possible to 
directly measure this self-coupling, it plays a role in how couplings run at higher energies and 
so indirectly constrains the possibilities for UV physics. 
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Figure 11: Constraints on Higgs-bulk mixing from bulk-nucleon bremsstrahlung in SN1987a. The 
left panel shows the constraints obtained assuming Tg„ = 20 MeV and the right panel shows the 
analogous constraints assuming TgN = 60 MeV. Couplings are evaluated at l/f — nih — 125 GeV. 



These constraints differ in the presence of Higgs-bulk mixing, because this mixing changes 
the relation between the Higgs mass and its self-coupling. Furthermore, brane-bulk interac- 
tions quite generally introduce a new source of running for brane-localized interactions. Both 
of these observations work to change the nature of the constraints on UV physics. 

There are generically two kinds of UV constraints that arise for Higgs couplings: 'vacuum 



stability' [56| and 'triviality' [57|. The first of these demands that the relevant quartic coupling 
of the Higgs potential not run to negative values. The second demands that any Landau poles 
(where the running couplings diverge) not arise at too low an energy. The depth of one's worry 
about these bounds depends fairly strongly on the depth of one's convictions as to how far 
this running can be trusted before some at-present-unknown UV physics intervenes. In the 
present instance triviality turns out to provide a fairly strong constraint on the brane-bulk 
coupling A2. 

To see how Higgs-bulk mixing changes things, we briefly restate the two sources of run- 
ning. For the Standard Model Higgs self-coupling, one-loop renormalization — including the 
Higgs, top quark and the gauge bosons in the loop — gives the following beta function for A 



(3.65) 



where yt = mt/v is the top-quark Yukawa coupling, and gi and g2 are (respectively) the U{1)y 
and SU{2)l gauge couplings, with the mass of the W and Z bosons given by = g'^v'^ 
and m| = | + gf) v'^ at tree-level, as usual. To these must be added the new Higgs-bulk 
contributions to the UV running of brane couplings, and writing /xd/d// = — /(f) d/df these 
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are given by 



^ dfi ) 27ra ' \ dfi) 27ra ' \ dfi J lua 





\ _ 5A2 















^d;ij, = w v^^J.^w y'dii)^ = ^- ('-^^^ 

where a derivation can be found in Appendix 

We've seen that C, and Acfr are two combinations of these couphngs appear quite frequently 
in Higgs observables. Eqs. ( p. 66 ) imply these couplings satisfy 



dC\ ^ ^CA2\ ]_ 

dnj g \TTaJ l + (A2/4a) 



/^^) =(^) . . ..„.2 > (3-67) 



dAeff A ^ n\ 1 - (A2/4a) 2 



and 

Notice that there is no requirement that A2 be small, so we need not expand the denominator 
in these expressions. 



Triviality 

An example of a solution to eqs. (|3.66|) is 



-9if^) = , , '^""'^ , (3.69) 

where we show the running relative to the Higgs mass scale, /i = m/j. This shows how 
strongly the running depends on the coupling A2. The triviality bound comes from de- 
manding that couplings like g remain within the perturbative regime throughout the energy 
ranges of interest. In particular, we require that the Landau pole (where g ^ 00 and so 
A2(n^/i) log(/u/m/i) = 27ra), not occur within this energy range. 

Notice that the Landau pole arises for scales < rrih whenever A2 < 0, but occurs for 
H > nih when A2 > 0. If we demand no such pole at energies below 1 TeV then we must 
therefore require A2 to lie in the range 

< A2 = — < 3 , (3.70) 
a 

as indicated in Fig. |l|. The upper limit becomes smaller if no Landau pole is allowed for 
energies above 1 TeV. This represents a significant constraint since none of our results required 
perturbing in A2, and so were not restricted a priori to small A2. 
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Vacuum stability 

The vacuum-stability bound demands that A remains positive as it is extrapolated into the 
UV (at least up to the point where any new UV physics intervenes to change how things 
run). In the Standard Model this provides the strongest constraint for light Higgs masses, for 
which A must be small. In this case, neglecting A in its RG equation gives 



(3.71) 



and the constraint arises because the dominant term is negative, eventually driving A negative. 

The presence of Higgs-bulk mixing can change this constraint in at least three ways, 
two of which act to weaken the vacuum stability constraint. It first does so by moving the 
starting point for A to more positive values. That is, the condition ruh = 125 GeV implies 
Aeff (m/i) = 0.1291 and so 

A(,nO = OA291 + |(|il^. (3.72) 

rather than simply \{mh) = 0.1291, as would have been the case for the Standard Model. 
Because the triviality bound requires A2 > we see that the presence of bulk couplings moves 
the initial condition, A(m/j), to more positive values. 

The second change is to the RG equations governing the running of A. Including both 
the contributions from the bulk and from Standard Model loops, we have 

l^f-) =^SM + ^, (3.73) 

where Psm represents the contribution of eq. ( |3.65| ). This shows that the bulk couplings 
always make dX/dfi more positive, and make it more difficult for A to become negative at 
higher energies. 

Finally, the third way Higgs-bulk mixing changes the logic of these bounds is by providing 
new UV physics, beyond which a naive extrapolation using the renormalization group need 
not apply. Perhaps the most dramatic way this might happen can be seen in the large- volume 
case, for which the KK scale is much smaller than the Higgs mass. In this kind of scenario 
the extra-dimensional Newton constant is much smaller than the 4D Planck mass, and all 
extrapolations must break down at the mass scale associated with the extra-dimensional 
Newton constant due to the intervention of the UV physics (perhaps string theory) that is 
required to unitarize gravitational interactions. 

Ultraviolet Fixed point for A 

The new, positive bulk contribution to the beta function also allows A to reach an ultraviolet 
fixed point. Demanding the vanishing of eq. ( p. 73 ) at some UV scale fi gives 



^ « [P + Ayf - yt - ^(3^1 + gf) + ^(Sffl + 2gh! + gf) ■ (3.74) 
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To approximate the numerical value of this expression we evaluate the gauge couplings, top 
Yukawa coupling and Higgs quartic coupling at the weak scale, which is justified by their 
weak running and the proximity of /i to the electroweak scale. (A scale much higher than the 
weak scale would exceed the extra-dimensional gravity scale, as discussed above.) We also 
approximate \{mh) = 0.1291 + 0{g'^{mh)) 0.1291 and so we find that g{^)/y/a 0.4 is 
required to obtain a UV fixed point. 

This condition can be run down to the Higgs mass scale using eq. ( p. 691 ), which gives 



gi^h) gin) 



a \ a 



27ra l°g(/^/"^^) 



(3.75) 



Since constraints on Higgs bulk mixing restrict g{mh) / y/a <^ g{n)/^/a ~ 0.4, the value 
of A2(?Ti/i) required to attain a UV fixed point is given approximately by the vanishing of 
eq. ( |3.75 ). Choosing A2(m/j) in this way also ensures a Landau pole at /i, and so we conclude 



that the UV fixed point for A and Landau pole are approximately coincident. 



4. Conclusion 



In summary, this paper examines the phenomenological implications of Higgs mixing with a 
bulk scalar field within an extra-dimensional brane-world scenario with the Standard Model 
localized on a brane. We focus in particular on the Higgs portal: the dimensionless couplings 
that can exist in such a scenario between a brane-localized Standard Model Higgs and a bulk 
scalar field if there are precisely two extra dimensions. We have a variety of motivations for 
studying this problem, with the main one being the requirement for such couplings in the 
recently discovered mechanism |22| for stabilizing two dimensions at naturally large values 
(through a manner similar to the Goldberger-Wise mechanism [^l| with one extra dimension). 

In a nutshell, we find that the Higgs portal causes the Higgs to mix with KK modes of 
the bulk scalar, generically leading to new channels for emitting missing energy in processes 
that can lead to Higgs emission. This can give observable signals at the LHC that strongly 
resemble the phenomenology of a Higgs with a branching ratio to an invisible decay channel. 
Unlike Higgs-curvature mixing, or bulk emission by a Higgs (through the trilinear hh(j) cou- 
pling), Higgs-bulk mixing through the extra-dimensional portal can give rise to an appreciable 
energy-loss rate in supernovae, leading to significant constraints on their couplings. As shown 
in Fig. [l|, although strong, these constraints need not preclude an observable signal for an 
invisible Higgs 'decay' channel at the LHC. 

In more detail, we find: 



• The strongest constraints come from nucleon bremsstrahlung in SN1987a: the phase 
space made available by the high temperature of the supernova, the appreciable Higgs 
coupling to nucleons, and the strongly coupled nature of nucleon interactions makes this 
a very strong bound, as we expect from experience with graviton emission in extra di- 
mensional models. However, the one pion exchange approximation we employ is known 
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to overstimate the emission rate, the temperature of the supernova core is uncertain 
and our calculation assumes the non-degenerate limit, which overestimates the rate 
at small temperatures |5^. Furthermore, the neglected effect of multiple scatterings 
will decrease the emission rate. Therefore we conservatively estimate bounds assuming 
TsN = 20 MeV, which gives the dominant bound g/\/a < 1.5 x 10"'^ (or ( < 5.6 x 10"'^) 
when A2 = 0. 

• LHC global fits: the LHC can indirectly probe an invisible Higgs. An additional, invisible 
width supresses Higgs signals, and too large a suppression would be in tension with the 
Standard-Model- like strength of signals currently being observed at the LHC. Global 
fits to the Tevatron and LHC Higgs data currently imply bounds i?inv < (0.3 — 0.64) 
the most conservative of which imposes g/\/a < 0.007 (or ( < 5 x 10^^), a bound 
subdominant to the SN1987a bound. If more LHC data were to begin favouring a 
universal suppression to Higgs signals, then this would suggest an invisible Higgs width, 
possibly from Higgs-bulk mixing. 

• Invisible final states at the LHC: the LHC can also directly search for invisible Higgs 
decays. The LHC will be most sensitive to the 2j + ^j, signal from vector boson fusion 
into an invisible Higgs. Current searches are only sensitive to invisible Higgs cross 
sections roughly equal to the Standard Model cross section. However, at design energy 
and 300 fb~^ of integrated luminosity, it was estimated that this channel would allow 
the 5a discovery of an invisible Higgs with Standard Model production cross section 
and an invisible branching ratio as small as -Binv ~ 0.25. With a more modest 20 fb~^ of 
luminosity at 7 and 8 TeV, the LHC should be able to rule out i?inv ^ 0.4 at 95% C.L. 
In both cases, there are regions of allowed parameter space that predict these branching 
ratios. 

• Additional cubic interaction: Higgs-bulk mixing also predicts a cubic Higgs-Higgs-bulk 
interaction ^gh'^4>{0). This interaction was studied in in the context of the h<p ^ 
'jjcj) final state before the discovery of the new 125 GeV Higgs-like resonance, and 
without A2 and Higgs-bulk mixing effects. It was concluded that the LHC can probe 
down to g = 0.18 with 100 fb~^ at 14 TeV. The possibility of large A2 and a h(p ^ bbcj) 
final state through this interaction make the 77 and bb signals worth investigating 
as a probe of Higgs-bulk physics. The observation of an invisible Higgs width consistent 
with an increased rate in one or both of these channels would provide strong evidence 
of Higgs-bulk mixing. 

• Beyond the large-volume limit: we present constraints and signals in the large-volume 
limit, but expect the the phenomenology to change a great deal if the volume were not 
large, so that ruh ^ rnj^i^. A single state in the diagonal KK tower would be identified 
as the 125 GeV resonance, and the spacing of nearby states would be governed by R~^. 
There is no reason to believe that the main signals of large R Higgs-bulk mixing — 
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invisible final states and suppressed Higgs production rates - would persist, and the 
astrophysical bounds would no longer apply to the much heavier states of this scenario. 
This scenario might resemble Higgs-radion mixing in Randall-Sundrum models, and we 
regard its exploration to be worth pursuing. 



Future colliders: although we refrained from discussing these in any detail, a future 
muon collider running at ^/s = rrih with 0.5 fb~^ of data and beam energy resolution of 
0.01% (0.003%) could directly probe the Higgs width to within 0.85 MeV (0.30 MeV) 



59 1 . A future linear collider with an integrated luminosity of 250 fb at 250 GeV 



would be able to constrain the Higgs invisible branching ratio to less than a few % |6C]. 
This would correspond to the B = 0.03 line in Fig. |l]. 



In short, mixing through the Higgs-bulk portal provides a particular example of what 
detailed Higgs studies might ultimately tell us about the nature of vacuum energetics. We 
hope that this is the path Nature chooses, and that the Higgs is shown to have exotic invisible 
properties, rather than persisting in having invisible exotic properties. 
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A. Renormalization of brane couplings 
Codimension-1 Regularization 

In this Appendix we compute the near-brane divergences and how they are renormalized. 



following the analysis of Ref. |19|. 



A more formal means of deriving the regulated boundary condition in ( p.l5| ) and the sub- 
sequent renormalized brane couplings ( p.l6|) is to model the codimension-2 brane as the dimen- 
sional reduction of an arbitrarily small codimension-1 brane. We introduce a codimension-1 
brane action at a small distance r = e in the extra dimensions, such that its dimensional 
reduction matches the codimension-2 brane action as e — )• 0. This might not seem like a 
palpable addition to the theory. However, since new brane has codimension 1, it resolves the 
singularities typically associated with codimension-2 sources, at the expense of an introducing 
arbitrarily short distance e. It is in this sense that the codimension-1 brane is analogous to 
introducting cutoff in typical quantum field theories. 
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Consider the following codimension-1 action 

S, = -f{e) J d^xd9(T,-{fi%fH^H + X,{H^Hf + {fi%)^'^, + ^Xl<^^,+g,H^H<P,^ , 

(A.l) 

where the factor /(e) comes from the determinant of the induced metric, = ^{x,e,9) is 
a function of x and 9 evaluated at r = e, and the rest of the Standard Model terms are 
neglected. Couplings with an e script are new couplings defined on the codimension-1 brane, 
in analogy with the quantities in (2^). The field H is not a function of the extra dimensional 
coordinates. 



We can match the codimension-1 action to the codimension-2 action in (2.4) by dimen- 
sional reduction. The KK modes of $e in the 6 direction acquire a mass n/e so all but the 
n = modes decouple for arbitrarily small e. The surviving n = modes have constant 
profiles in the coordinate 6 that we can choose to be unity, giving 

= -2naf{e) j d^x (t, - (lil^fH^H + K{H^ Hf + + + g.H^H^^ , 

(A.2) 

where we have identified $b = <^(x,0) as the zero mode of To reproduce the action ( |2.4D 
requires 

27ra/(e)c, = c, (A.3) 

where is any of the couplings in the codimension-1 action (STl) and c is the associated 
coupling in the codimension-2 action ( [2.41 ). 

Having matched the codimension-1 action to the original theory, we use it to derive a 
boundary condition for The variation of the bulk action (2T) gives a contribution on the 
boundary at r = e 

6SsD- J d^x de 6<^, i-fdr^X^, , (A.4) 

which must vanish when combined with the variation of the codimension-1 action (|A.lD , 
giving the boundary condition 

-{fdr^)r=e + /(e) ((^1)' + A^2^, + g^H^ = . (A.5) 

This can be rewritten in terms of the couplings in the codimension-2 action 

-27ra ifdr^)^^^ + /i| + Aa^e + gH^H = , (A.6) 

which is a regularized version of eq. ( |2.9| ). Varying (|A.1[) with respect to H yields the 
regularized version of ( 2.11| ), which can be used to derive ( 2.15| ). 

Floating Brane Renormalization 

As in typical field theories, the cutoff e is not a physical quantity, and as e the boundary 
condition still diverge. To remedy this, we consider a theory that is sourced by a "fioating" 
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codimension-1 brane at r = f, that can be dimensionally reduced to a codimension-2 brane. 
The properties of the floating brane are fixed by demanding that the solution to the field 
equations for r > f matches the regularized solution. The dependence on e will be eliminated 
when the bare couplings c are traded for the coupling constants of the floating brane c. It is 
in this sense that f can be thought of as a subtraction scale, and c can be thought of as the 
renormalized couplings. 

After dimensional reduction to codimension-2, we assume the floating brane has the form 



Sb 



J d^x (f - fllH^H + \{H^Hf + Jll^r + \\2^l + gH^H^r^ , (A.7) 



where = ^{x,r) and the theta dependence has been integrated out. Varying the action 
gives two boundary field equations 

-27rQ ifdr^l^, + gH^H + Aa^, + 7l| = 0; H = i - g^,) , (A.8) 

2a 



in direct analogy with the boundary field equations in §2.2. These equations should be read 
as fixing the floating brane couplings, since we have demanded that $ and H solve the 
regularized field equations and so their functional form is already fixed. For a given value 
of f, the floating brane couplings will have to be chosen appropriately, and will change with 
a change in f. However, since f is arbitrary, this change in the couplings cannot have any 
effect on physical quantities, such as H^H. For example, changes in ( |A.^ ) under a change in 
f should vanish, giving 

H^Hdrg + $(f)5f A2 + AaSf $(f) + drjil = 

2H^Hdf\ - dfUl + ^{f)dfg + gdr^{f) = , (A.9) 

where we have used the fact that df{fdr^)r=r = {drfdr^)r=f = by the bulk equation of 
motion (|2.6D , and df{H^H) = because it is a physical quantity that should not depend on 
f. We premultipy both equations by /(f) to facilitate the use of the relation 

f{r)drHr) = {fdr^)r=f = ^ (a^I + ^2$(f) + gH^H^ . (A.IO) 

Substituting this into ( |A.S| ) and equating powers of <I> and H'^H yields the following RG 
equations 

r. _2 ^2pq> r. _ g^2 r, T ^2 

l>f^ = a,4 = g; 9,f = -A, ,A.n) 

where we have used f{r)df = dp. The solutions can be found in ( p. 16 ). Although not 



explicitly derived in this section, for completness the RG equation for T is listed here and 
can be derived from ( |B.5D in the following appendix. 
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B. Probe-brane energetics 



In this section we evaluate the energy density to identify when the solution with nonzero 
H is energetically preferable to the solution H = 0. For illustrative purposes (following the 
discussion of |jl^) we do so here ignoring the energetics of the gravitational back-reaction 
of the branes. This is known not to be a good approximation in general for codimension- 
2 objects since in many cases it is the classical bulk back-reaction that cancels the brane 
tensions to allow the classical brane geometries to be flat |24, 30, (Indeed, it is this 



property — together with bulk supersymmetry — that underlies their use as a starting point 



for seeking solutions to the cosmological constant problem 27, p2[.) 
The regularized energy density for H and $ is 



n = ns + Y.Ub, (B.i) 



where the bulk contribution is 



ttR 

2 



nn = 7ra / dr/(r)(9,$)^ (B.2) 



and the 'Higgs brane' gives 



U = T-fil H^H + X{H^Hf + iil^, + ^^l+g H^H . (B.3) 



Evaluating using the bulk and brane solutions, eqs. ( p. 19 ) and ( 2.20| ), and renormalizing 
as before gives 



47ra J V 4vrQ! J i , ^^jr) p/ 4A(f) 

which suggests we renormalize the brane tension using 

T{f) = T+i^\ ^'"^ , (B.5) 

SO that the vacuum energy density coming from the bulk and the Higgs brane is 

n. + U = - (%[^) F(., .R) + m - §1 . (B.6) 

Eq. ( |B.6| ) allows two simple conclusions to be drawn. First, since the result with H = 
only omits the last term, this shows that nonzero H is a preferred vacuum to the extent that 
Ji\{r) > 0. (Boundedness oiTi from below precludes A from being negative.) Second, the sign 
of the p.%gg term depends on the sign of F(r,7ri?), which in turn depends on the boundary 
condition at ttR . We can imagine imposing a general boundary condition 

27raf<^'{7rR) + \3<^{irR) = , (B.7) 
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which fixes f as foUows 

Fif,7rR) = ^. (B.8) 

This illustrates that the bulk field contribution to the vacuum energy density can be positive 
or negative depending on the sign of A3, so that the preferred vacuum will depend on the 
details of the faraway brane. We assume that the faraway brane is such that the general 
solution we found in the main text is preferred. 



C. Schwinger-Dyson equation 

In this Appendix we derive the relation of eq. ( [3.24 ) 



[Amp{h<l)*{0))k]* = Ampim h*)k 



ig{r)v {hh*)k 



l-i(A2(f)/4a) 



(C.l) 



f2 = _l/fc2 



relating the amputated mixed h — (p propagator to the h — h autocorrelation, and the relation 
of eq. ( pSD 



{hh*)k 



l+[iX2 +M2Z)^]Df(0,0) ' 



(C.2) 



that gives the dressed two point function h, and largely controls the phenomenology of Higgs- 
bulk mixing. 

Our goal is to compute relations amongst the four correlation functions of interest, given 

by 

{hh*)k := Ghh{k) , {h<P*{y))k := Gh^Mv) , 

{Hy)h*)k ■■= G^h(.k;y) and {4>(,y) (t>* (,y'))k ■= G^^{k;y,y') , (C.3) 

where Gh(j,{k;y) = G*^f^{k;y). In these expressions denotes the spatial coordinates in 
the two extra dimensions while k^ is the Fourier transform variable in the four on-brane 
directions. 

The most direct way to obtain the desired relations is to express the Higgs-bulk interac- 
tions as delta-function localized terms in the lagrangian density, following arguments made in 
the appendix of ref. [0].^ The starting point is the field equations for linearized fluctuations 



g2(n4 + n2) 



\2<P + gvh 



5\y) = 



□4/1 - 2At;^ h- gv (/)(0) = . 



(C.4) 



®A disadvantage of the delta-function formulation is the requirement to deal with expressions like f{x) 5{x), 
with f(x) CX3 as a; — > 0. This requires a more careful treatment of regularization and renormalization, along 
the lines of the codimension-one formulation used in the main text, but in the present instance leads to the 
same conclusions. 
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which imply the fohowing equations for the propagators 

V^(-fe2 + U2)G^^{k; y, y') - [x^ G^^{k; 0, y') + gv Gh^k; y')] 6^{y) = iS\y - y') 

V^(-A;2 + □2)G^h(fe; y) - [\2 G^hik; 0) + gv Ghh{k)] 5\y) = 

{k" + 2\v^)Ghh{k) + gv G^hik; 0) = -i 
(k^ + 2Xv^)Gh4k; y) + gv G^^k; 0, y) = . (C.5) 

By contrast, the unperturbed propagators in the absence of Higgs-bulk couphngs satisfy 



(C.6) 



g2{-e + n2)Dt{y,y') = i5\y-y') 
{k'^ + 2\v^)Dl = -i. 



We use the first of eqs. ( C.6 ) to solve the second of eqs. 



leading to 

G^h{k;y) = -ij d^y'Dl{y,y')\\2G^h{k;^)+gvGhh{k)Y{y') 
= -iDtiy, 0) [A2 G^hik; 0) + gv Ghh{k)] , 
and this, when specialized to y = 0, in turn implies 

G^h{k; 0) = -iDtiO, 0) [A2 G^h{k; 0) + gv Ghh{k) 
which may be solved to give 



(C.7) 



(C.8) 



G^h{k;Q) 



gv Ghhjk) 
1 + 1X20^0,0) 



DIM. 



(C.9) 



The overall factor of Dt (0,0) is removed when the external (/)-line is amputated, and for 



the denominator we use the continuum result, eq. ( p.36| ), to evaluate 1)^(0,0), 



Ana 



log(— /c^e^) — in 



(C.IO) 



and renormalize the divergence into the brane couplings, g and A2, using eqs. ( 2.16 ). Eq. ( p.l| ) 
then follows by choosing the renormalization point so that /c^f^ = 1 and the logarithms vanish. 



We can also use the second of eqs. ( p.6|) to solve the third of eqs. (|C.5D , giving 

Ghh{k) + igvG^h{k;0)D^ = D^ (C.ll) 



To solve for Ghh{k) we substitute eq. (|C.9D into the above expression, which can be rearranged 
to give eq. (|C.2|) as desired. 
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D. Toy model: unperturbed modes 



In this appendix we explicitly take the continuum limit of ( 2.33| ) to arrive at eq. ( pBI) . This 
is accomplished in a toy model in which the extra dimensions are a flat disc: /(r) = r for 
< r < ttR. We can explicitly solve the wavefunctions on this background, which allows for 
a straightforward move to the large R limit, although the results are true for all R. 

Using eq. ( 2.29| ) in eq. ( p. 28 ) on the disc geometry gives the field equation for the n = 
wavefunctions 



M^l + -dr{rdr) 



with the following boundary conditions 

irdrPoi)r=o,nR = 0; 

and normalization conditions 



27ra J drrP^iPoi> = Sii> . 


The properly normalized solutions and eigenvalue conditions read 

1 



(D.l) 



(D.2) 



(D.3) 



Poiir) 



V^^3^ \jQ{MonTR) 



with Ji{7rRMoi) = 0, 



(D.4) 



where Jq is the zeroth Bessel function of the first kind. The brane-to-brane propagator is 
given by 

-i \ MeMoi) 



^ , A;2 + m2 -iej J^{7rRMoi)7r^aR^ ' 
since Jo(0) = 1. Using the fact that Jo{x) — for large x gives 

Jo{eMoi)Moi 



(D.5) 



2TTaR 



(D.6) 



Sums over closely spaced modes in d dimensions can be replaced by integrals as follows 



(D.7) 



where, in this case, the sum is over the radial index, so the conversion is one-dimensional. 
Using il*^ = 27ri? for the diameter of the disc is gives 



Dt{e,0) 



lira 



dq 



qJojeM) 

^2 _|_ g2 _ 



(D.8) 
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The integral can be computed 



where Kq is the zeroth modified bessel function. For small arguments Kq(x) — )• — log(x/2)+7 
so the divergent part of the brane-to-brane propagator reads 



^f(e,0) = ^log(feV), (D.IO) 



in agreement with eq. ( 2.36| ). 



E. Beyond Sturm Liouville 

In this Appendix we describe how the Sturm-Liouville orthogonality conditions generalize to 
the case of interest in the main text, for which the boundary conditions differ for different 
modes. 

For the present purposes the eigenvalue condition for the mode functions (,nix) has the 
general form 

dx [p{x)dxCn] - q{x)in + Kw{x)Cn = , (E.l) 

in an interval xq < x < xi, with p,q,w known real functions and A„ the corresponding 
eigenvalue. The unusual part relative to Sturm-Liouville problems of childhood days is that 
they satisfy n-dependent boundary conditions at the edges of the domain of interest: 

JbiXn - Kh)pd^in + (A„ - Lb)eJ = , (E.2) 

- X=Xf, 

where J^, Kf^, Lf, are again known coefficients. These boundary conditions ruin the orthonor- 
mality of the mode functions under the usual inner product, 

XI 

dx w{x) CmCn / Smn , (E.3) 

XQ 

which in turn ruins the diagonalization of the 4D action once decomposed in terms of these 
modes. 

To identify how the inner product must generalize in order to maintain orthogonality 
with the new boundary conditions we follow standard steps. First multiply ec[. ( [E.l[) by 
then subtract the complex conjugate of the same equation with (m -f-)- n) and integrate the 
result over x. This yields 

XI 

{Xm - A„) / dx W{x) CmCa = \p {Cm dxU " Cn d.Cj] , (E.4) 



XQ 



XO 
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which would vanish for the usual Stum-Liouville boundary conditions. However, with the 
n-dependent boundary conditions of the form ([E.2D we instead have 



Xl 

v„)/< 



(Am - A„) / dx w{x) il^^n = {K 



b 



Lb-Kt 
Jb 



Cmi^b)^n{xb) 



(E.5) 

What allows us to devise an inner product with respect to which the modes are automat- 
ically orthogonal is the property that the n-dependence of the boundary conditions is linear 



in Xji, since this ensures both sides of eq. (E.5) depend on n through their common factor of 
(Am — A„). This suggests defining the following inner product 



Xl 



Cmixb)^nixb) 
(Am — Kf,){Xn — Kb) ' 



xo 



since eq. (E.5) then shows that the boundary conditions imply modes with different eigenval- 
ues are automatically orthogonal, and so a basis of eigenmodes can be chosen to be orthonor- 

mal: {^mi^n) = ^mn- 

In the dimensional-reduction problem the constants Jb,Kb,Lf, are read from the brane 
action, and so are the quantities that appear in the quadratic lagrangian once bulk fields are 
decomposed in terms of these mode functions. This ensures that the action diagonalizes as 
it would have done for a standard KK decomposition without endpoints. 

For example, for the zero modes in the brane bulk mixing scenario we send x — )• r and 
n, m — )• s, t and use p{r) = 27ra/(r), q{r) = and w{r) = 27ra/(r). We replace the eigenvalues 
with the KK masses A^ = Mg^. There is only one brane at ro = and (neglecting subscripts) 
it gives J = -I/A2, K = 2v'^X and L = 2v^X + {gvf/X2 so that the inner product reads 



{rs,rt)=27ra / drfVtVt + 



(M2 - 2Xv^){Ml - 2Xv^) 



(E.6) 



and the orthonormality relationship ( ^.50 ) in the text follows 



Diagonalization of the quadratic action 

We now show that in the case of interest in the main text, this modified inner product is just 
what is required to diagonalize the quadratic action, including the Higgs-bulk mixing terms. 
For simplicity, we only include the n = modes, but the extension to any n 7^ level of the 
KK tower follows readily. We still use s = {n, i} with the understanding that n = 0. 
In terms of KK modes the the bulk action ( p. 22] ) reads 



Sb = -2vrQ / d^x I drY.UK'Pt) 

s,t 

-27rQ J d^x J dr^ if drV*drV, 



(E.7) 
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where terms have been organized into their r-dependent parts, which are in round brackets, 
and their x-dependent parts, which are in square brackets. They have also been written on 
separate hues for organizational purposes. Integrating the second term by parts gives 

^1 



Sb = -2vra / d^x j drY.if'Pl'Pt) 

•' •' s,t 

+2-na I d^x j dr^^ {Vt drfdr V, 

s,t 

+2™ / d^xY,imdrv:),^o 

s,t 



where the term on the bottom line is a boundary term, and it is assumed that the other 
boundary term for the faraway brane vanishes by the boundary conditions, or is cancelled by 
the faraway brane's action. The bulk equation of motion (2.45) allows the second line to be 
combined with the first as follows 



Sb = -27Ta I d^x I drY^UKT^t) 

s,t 

+2™ / d''xY{fVtdrV:\^^ 



1 



(E.9) 



The integration over the radial coordinate can be completed using the orthonormality rela- 
tionship ( 2.501 ) so that the bulk action contributes three terms to the dimensionally reduced 
Lagrangian that will be called ^1,2,3 



JO2 



E 



{gv)^Vs{0)Vt{0) 
(M| - 2At;2)(M2 - 2Aw2) 



E 



C3 = 2TiaY,{mdrVs),-_ 



s.t 



(E.IO) 



The first line is a canonically normalized KK tower of scalar fields with masses Mg, which 
is the desired final result. The second and third term are cancelled by terms in the brane 
action as will be shown explicitly. For example, writing the h kinetic and mass term in terms 
of eigenstates <ps and then combining them with C2 gives 

{gv)^Vs{0)Vt{0) 



1 



2l2 



C2 + Ch = C2- -df.hd^'h - Xv'h 



s.t 



(M2 - 2A^;2) 



1 



while and the brane mass term for (j) give 



s.t 



2 'Ps^t 



(E.ll) 



(E.12) 
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which is identical to ( |E.ll ) once the boundary condition ( 2.49 ) is employed. Finally the 
mixing term gives 



C 



(E.13) 



so that £2 + + £0 + £/i + Ch(f) = and the dimensionally reduced theory is the KK tower 
of massive scalar fields found in Ci. 



F. Solving for the KK mode functions 

In this appendix we solve solve the n = perturbed wavefunctions Vqi for /(r) = r, and 
< r < ttR with Dirichlet boundary conditions at r = ttR. In this geometry the general 
solution can be written in terms of Bessel functions 

Voeir) = Ne [|lo(rMo^) + DeMrMoi)] , (F.l) 

where Ni are normalization constants, are integration constants and the factor of 7r/2 is 
chosen for convenience. It is straightforward to impose Dirichlet BCs at r = ttR, which imply 

7t Yo{7tRMo,) 

^ 2 JoiTTRMoi) ■ ^ ' ' 

Imposing the UV boundary condition, on the other hand, is more complicated. This is because 
of the UV divergences we expect in this theory. Near the origin, the relevant Bessel functions 
behave like 

Yo{x) « -[ln(x/2)+7] 

TT 

9 1 

yi(x) « --- 

vr X 

MO) = 1 

Ji(0) = 0. (F.3) 

So, as in the vacuum solutions, the boundary condition near the brane diverges, and must be 
regulated and renormalized. We cut off the boundary condition at r = e and find 

De = ^-log{eMoi/2)-j with /3, = A2 + -^^^^^ . (F.4) 



We can rewrite this boundary condition in terms of the renormalized quantities of ( p. 16 ), 
rendering it finite and cutoff-independent 

A = |^-log(r-Mo,/2)-7 with -p,{f) = \2{f) + , ./^"^f.,-, • (F.5) 
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Equating the two expressions for Di yields an eigenvalue equation for the M| masses 

TrYoiirRMoi) _ 27ra 



2 JoiTrRMoe) ^^{f) 



log(fMo£/2) - 7 , 



(F.6) 



which is, unfortunately, quite difficult to solve. 

The normalization condition for the perturbed wavefunctions reads 



TVR 



(M2, - 2v^\f 



1. 



(F.7) 



We break this calculation into parts. First, we calculate the integral 

■kR ttR 



I 



rdr Vnp (r) = No I rdr 



l2 



Using the identity 



d_ 

dx 



xZq{x) , 



(F. 



(F.9) 



for any function Zq that satisfies Bessel's equation, we can write 



irRMn 



2 



(^^YoinRMoe) + AJo(vri?Mo,))' + (|yi(^i?Mo,) + DiJi{^RM^ 



2M2 ' 



where the second term follows from taking the e — )• limit. Using the boundary condition 
(13) gives 



, 1 7r2/?2 
Ip = -- 



IT 



-Yi{7rRMoe)Jo{7rRMo^) - -yo(vri?Mo^) Ji(7riiM, 



2J^{7TRMoi) \2 
Bessel functions obey the following identity 



^Yi{TTRMi)Jo{TTRM,) - ^Yo{TTRMi)Ji{TTRM^ 



1 



1 



2M2 



■nRMi ' 



so that 



. (F.IO) 



(F.ll) 



(F.12) 



2M2, \J^{t:RMm) 

Now we move to the second term in (|F.7D , which is equal to Bj. Using the r = boundary 
condition allows us to write 

gvVQi{Q) _ 2iTagv [rdrVpel^^Q 2-nagvNe, 



M2 - 2v^X (M2 - 2v^X)X2 + {gvY {M^„ - 2v^\)\2 + {gvY ' 



(F.13) 
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which can be inserted into (F\7) to give the following equation 





(F.14) 



[A2 {Ml,-2v^\) + {gvYf 



From the above normalization we find the mixing coefficients 



.-2 



1 + 




- 1 



) 



[A2 {M^,-2v^\) + {gv ff 



(F.15) 



2M2, \J^{^RMq,) 



2'Ka{gv)'^ 



Note that the mixing coefficients vanish as 5 — t- or A2 — 00 unless = 2v'^X 
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